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Abstract 


Current  advances  in  microfabrication  allow  for  implementation  of  miniaturized  electro-mechanical 
systems  that  fit  on  a  single  integrated  circuit  die.  Microeiectro-mechanical  systems  (MEMS)  have  many 
practical  benefits,  including  integration  with  control  and  signal  processing  electronics,  batch-fabrication, 
and  low  cost.  Moreover,  the  small  individual  mass  of  the  micromachined  devices  leads  to  superior 
ruggedness  and  fast  system  response  time,  making  MEMS  ideal  for  a  variety  of  military  and  commercial 
applications. 

One  method  for  MEMS  manufacturing  is  surface  micromachining.  The  challenge  with  this  method 
is  its  inability  to  produce  structures  with  high  aspect  ratio  due  to  its  planar  fabrication.  The  surface 
micromachined  structures,  however,  can  be  assembled  out  of  the  plane  of  the  substrate  to  achieve  a  3- 
dimensional  function.  A  common  solution  for  assembly  of  highly  3-dimensional  MEMS  is  the  fabrication 
of  hinged  components  that  can  be  lifted  or  ”popped-up”  into  assembled  structures.  The  hinged  structures 
can  then  be  combined  with  microactuators  for  positioning  and  operating  MEMS.  Hinged  devices  are 
commonly  assembled  manually,  a  time-consuming  and  delicate  process  which  leads  to  low  yiekl  and 
inability  to  produce  devices  in  large  quantities.  The  research  objectives  of  our  work  focused  on  elimination 
of  such  fabrication  limitations  by  using  two  novel  self-assembly  techniques  for  MEMS:  assembly  using 
integrated  micromechanisms  (referred  to  as  MEMS  Self-Assembly),  and  assembly  using  a  solder  reflow 
process  (referred  to  as  Solder  Self-Assembly),  The  self-assembly  of  MEMS  eliminates  the  need  for 
manual  assembly  or  adjustment,  thus  making  batch  fabrication  feasible.  In  addition  to  manufacturability, 
the  self-assembly  of  MEMS  may  be  particularly  important  for  applications  which  require  reliable 
deployment  and  remote  assembly  of  delicate  structures  in  the  operating  environment  or  readjustment  of 
components  to  align  a  system  for  optimum  performance.  Specifically,  this  work  was  divided  into  eight 
subsections: 


•  The  study  and  development  of  self-assembled  comer  cube  reflectors. 

•  The  study  and  development  of  self-assembled  3D  switches. 

•  The  study  and  development  of  3D  micro-robot  legs. 

•  Design  for  MEMS  mechanism  driven  self-assembly. 

•  Design  for  surface  tension  driven  self-assembly. 

•  The  study  and  development  of  hybrid  assembly  structures. 

•  The  study  and  development  of  advanced  stmctures. 

•  The  quality  and  reliability  of  self-assembled  MEMS. 
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(Wx  W,  W=300  Jim)  and  partially  covering  gold  pad  (Lx  L):  (a)  L=180  jim,  (b)  L==250  jim . 42 

Figure  T8-4-6:  Experiment  result  shows:  (a)  the  L=250  jim  strucmre  bifurcates  earlier  than  the  L=180  Jim  structure,  (b) 

the  L=250  jim  structure  bifurcates  earlier  than  the  L=300  jim  structure . 43 

Figure  T8-4-7:  Prediction  results  of  the  cooling  process  from  100  "C  to  room  temperature  for  square  poly  silicon  plates 
(300  Jim  X  300  Jim)  with  square  gold  pads  (LxL).  (a)  L=100,  150,  175,  200,  and  225  |im;  (b)  L=225,  250, 275,  and  300 

Rm. . 43 

Figure  T8-5-1:  Curvature  variation  during  the  thermal  cycle  process  between  190  °C  and  room  temperamre  for  the 

center  supported  gol d/poly  l+poly2  200  Jimx  200  jim  square  plate . 44 

Figure  T8-6-1:  Measured  curvature  development  during  cooling  and  a  subsequent  stress  relaxation  process  at 
elevated  temperature  of  a  300  jim  x  300  jim  gold  (0.5  jim  thick)  /  polysilicon  (1.5  jim  thick)  plate  structure:  (a) 
Curvature  development  during  cooling  from  190  “C  to  120  ®C;  (b)  Curvature  development  due  to  the  stress  relaxation 

when  the  structure  was  held  at  a  constant  temperature  of  120  ®C . 45 

Figure  T8-6-2:  Measured  and  predicted  curvature  development  during  isothermal  stress  relaxation  process  at  elevated 
temperature  of  120  “C.  (a)  200  Jim  x  200  jim  gold  /  polysilicon  plate  structure  with  different  thickness  of  polysilicon; 

(b)  300  Jim  X  300  jim  gold/polysilicon  plate  structure . 45 

Figure  T8-7-1:  (a)  Multilayer  structures  after  indium  deposition;  (b)  schematic  diagram  of  the  plate  and  substrate,  red 
dots  indicate  the  locations  of  the  dot  test,  red  lines  indicate  the  locations  of  the  line  scan  test . 46 


V 


Figure  T8-7-2:  (a)  SEM  picture  of  the  line  scan  area,  (b)  Line  scan  result  which  shows  Si,  Au,  and  In  element 

component  distribution  along  the  scan  line . 46 

Figure  T8-7-3:  (a)  SEM  picture  of  the  dot  test  position,  (b)  Dot  test  results  which  show  element  component  at  the 

point . 47 

Figure  T8-7-4:  Dot  map  results  which  show  element  distribution  at  the  area . 47 

Figure  T8-8-1:  (a)  5  mil  solder  balls  were  put  on  top  of  circular  structures  and  rcflowed  ;  (b)  SEM  picture  of  room 

temperature  solder  shape  after  reflow . 48 

Figure  T8-8-2;  Measured  and  predicted  deformation  of  polysilicon  plate  with  solder  in  the  center,  (a)  Measured;  (b) 
Predicted . 48 


VI 


1.  Summary  of  Research  Objectives 


Current  advances  in  microfabrication  allow  for  implementation  of  miniaturized  electro-mechanical  systems 
that  fit  on  a  single  integrated  circuit  die.  Micro -electro -mechanical  systems  (MEMS)  have  many  practical  benefits, 
including  integration  with  control  and  signal  processing  electronics,  batch-fabrication,  and  low  cost.  Moreover,  the 
small  individual  mass  of  the  micromachined  devices  leads  to  superior  ruggedness  and  fast  system  response  time, 
making  MEMS  ideal  for  a  variety  of  military  and  commercial  applications. 

One  method  for  MEMS  manufacturing  is  surface  micromachining.  The  challenge  with  this  method  is  its 
inability  to  produce  structures  with  high  aspect  ratio  due  to  its  planar  fabrication.  The  surface  micromachined 
structures,  however,  can  be  assembled  out  of  the  plane  of  the  substrate  to  achieve  a  3-dimensional  function.  A 
common  solution  for  assembly  of  highly  3-dimensional  MEMS  is  the  fabrication  of  hinged  components  that  can  be 
lifted  or  "popped-up"  into  assembled  structures.  The  hinged  structures  can  then  be  combined  with  microactuators 
for  positioning  and  operating  MEMS.  Hinged  devices  are  commonly  assembled  manually,  a  time-consuming  and 
delicate  process  which  feads  to  low  yield  and  inability  to  produce  devices  in  large  quantities.  The  research 
objectives  of  our  work  focused  on  elimination  of  such  fabrication  limitations  by  using  two  novel  self-assembly 
techniques  for  MEMS:  assembly  using  integrated  micromechanisms  (referred  to  as  MEMS  Self-Assembly),  and 
assembly  using  a  solder  reflow  process  (referred  to  as  Solder  Self-Assembly).  The  self-assembly  of  MEMS  eliminates 
the  need  for  manual  assembly  or  adjustment,  thus  making  batch  fabrication  feasible.  In  addition  to  manufacturability, 
the  self-assembly  of  MEMS  may  be  particularly  important  for  applications  which  require  reliable  deployment  and 
remote  assembly  of  delicate  structures  in  the  operating  environment  or  readjustment  of  components  to  align  a  system 
for  optimum  performance. 
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2. 


Technical  Status  Report 


2. 1  Status  of  Effort 


For  reference,  the  schedule  of  tasks  as  set  prior  to  the  start  of  this  contract  is  included  below.  The  shaded 
time  period  blocks  indicate  completed  tasks.  In  the  remainder  of  this  report  Section  2.2  summarizes  the  completed 
work,  and  Section  3  lists  publications/consultations  related  to  this  research. 


Task  1:  3-D  Corner  Cube  Reflector 

Apr-98 

Apr-99 

Apr-00 

Apr-01 

a.  Design  actuator/micromechanism  for  MEMS  self-assembly 

b.  Demonstrate  a  MEMS  self-assembled  reflector 

c.  Develop  fluxless  solder  technology 

d.  Develop  a  quick  solder  deposition  technique  using  spheres 

e.  Demonstrate  a  solder  self-assembled  reflector 

f  Demonstrate  a  reflector  assembled  using  MEMS  and  solder 

g.  Characterize  the  performance  of  the  reflector 

Task  2:  Array  Of  3-D  Switches 

Apr-98 

Apr-99 

Apr-00 

Apr-01 

a.  Demonstrate  a  solder  seif-assembled  microswitch 

b.  Demonstrate  an  array  of  solder  self-assembled  microswitches 

c.  Characterize  the  performance  of  the  microswitches 

Task  3:  Array  of  3-D  Microrobot  Legs 

Apr-98 

Apr-99 

Apr-00 

Apr-01 

a.  Demonstrate  a  MEMS  self-assembled  microrobot 

b.  Control  gas  flow  effect  of  the  fluxless  soldering 

c.  Improve  solder  deposition  technique  using  an  evaporator 

d.  Demonstrate  a  solder  self- assembled  microrobot 

e.  Demonstrate  a  microrobot  assembled  using  solder  &  MEMS 

f.  Characterize  the  performance  of  the  microrobots  assembled 

g.  Determine  the  limitations  of  MEMS  assembled  microrobots 

h.  Improve  deposition  techniques  for  2-um  solder 

r 

i.  Demonstrate  a  soldered  microrobot  using  <10  um  solder 

WM 

j.  Explore  the  use  of  non-solder  materials  for  microrobots 

k.  Demonstrate  a  microrobot  assembled  using  MEMS  &  solder 

{removed  during  project  kick-off  mcolinc) 

■ 

1.  Characterize  the  performance  of  the  microrobots  assembled 

(removed  during  project  kick-off  rncciing) 

■ 

Task  4:  Design  For  MEMS  Self-Assembly 

Apr-98 

1  Apr-99 

Apr-00 

Apr-01 

a.  Develop  a  model  of  thermal  microactuators  for  mechanisms 

b.  Develop  1 -hinge  mechanisms  using  thermal  microactuators 

c.  Experimentally  determine  the  actuator  model’s  limitations 

d.  Experimentally  verify  the  1 -hinge  mechanisms 
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Develop  5-hmge  mechanisms  using  thermal  microactuators 
(removed  during  project  kick-off  meeting) 


Experimentally  verify  the  5 -hinge  mechanisms 
(removed  during  project  kick-off  meeting) 


Task  5:  Design  For  Solder  Self-Assembly 


Apr-99 


Apr-00  Apr-01 


Determine  the  limitations  in  scales  100  to  25  um. 


).  Develop  a  model  to  design  1 -hinge  self-assembly 


c.  Experimentally  verify  the  1 -hinge  model’s  accuracy 


Investigate  the  use  of  non-solder  materials 


e.  Determine  the  limitations  in  scales  <  25  um. 


f  Expand  the  year- 1  model  for  different  solder  materials 


g.  Develop  a  model  to  design  5-hinge  assembly 


1.  Conclude  the  study  of  non-solder  materials 


i.  Determine  the  limitations  of  controlled  alignment  sequence 


J.  Experimentally  verify  the  5-hinge  model 


k.  Develop  and  verify  a  model  for  controlled  alignment  sequence 


Task  6:  Hybrid  Technology 


Apr-00  Apr-01 


Determine  the  effect  of  soldering  on  MEMS  self-assembly 


Define  integrated  design  procedure  for  the  hybrid  technology 


Integrate  solder  and  MEMS  design  models 


Demonstrate  hybrid  technology  using  innovative  MEMS 


-Task  7:  Futuristic  MEMS 


Apr-99 


Explore  new  MEMS  soldered  with  over  1 2  hinges  (House) 


Explore  new  MEMS,  e.g.  doors,  using  micromechanisms 


Explore  new  MEMS,  e.g.  a  2-bedroom  house 


Explore  new  MEMS,  e.g.  a  3-bedroom  house 


Task  8:  Quality  and  Reliability 


Develop  3-D  deformation  models  for  micromirrors 


Test  micromirrors  to  validate  the  deformation  models 


Test  hinged  structures’  force-deflection  relationship 


Integrate  the  relationship  into  the  MEMS-  and  solder  models 


Identify  single-hinge  structures’  failure  modes 


Identify  multiple-hinge  structures’  failure  mode 


Quantify  mechanical  failure  on  single-hinge  structures 


Develop  reliability  models  for  single-hinged  structures 


Couple  the  deformation  model  with  solder  effects 


Quantify  mechanical  failure  on  multiple-hinge  structures 


Develop  reliability  models  for  multiple-hinged  structures 
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2.2  Research  Summary 

Task  1:  3-D  Corner  Cube  Reflector  (CCR) 


Thermal  actuators  and  micromechanisms  have  been  used  in  the  MEMS  self-assembly  of  a  CCR  [1].  Figure 
Tl-1  shows  scanning  electron  micrographs  (SEMs)  of  an  assembled  CCR  (a),  a  vertically  deflecting  thermal  actuator 
(b),  and  a  horizontally  deflecting  thermal  actuator  array  of  micromotor  (c). 


(a)  (b)  (c) 

Figure  Tl-1:  Assembled  CCR  (a),  and  associated  thermal  actuators  and  micromechanisms  (b,c). 

With  reference  to  Figure  Tl-1:  the  assembled  mirrors  shown  in  (a)  are  initially  fabricated  laying  on  the 
substrate;  the  vertically  deflecting  actuator  shown  in  (b)  gives  the  mirror  an  initial  lift;  and  the  micromotor  shown  in 
(c),  but  out  of  view  in  (a),  pushes  the  mirror  the  rest  of  the  way  up  until  it  is  locked  into  place. 

We  have  developed  a  fluxless  solder  technology  and  a  quick  solder  deposition  using  solder  spheres.  Solder 
spheres  can  be  placed  across  gold  plated  pads  attached  to  microstructures.  The  solder  and  microstructures  are 
heated  to  the  melting  point  of  the  solder  while  in  an  atmosphere  of  nitrogen  gas  and  formic  acid  vapor  [2].  The 
nitrogen  gas  and  formic  acid  are  used  to  inhibit  and  remove  metal  oxidation,  respectively.  This  technique  was  used 
to  assemble  the  CCRs  shown  in  Figure  Tl-2.  The  diameter  of  the  63  %  Sn  /  37  %  Pb  solder  spheres  was  8  mil.  A 
vertically  deflecting  thermal  actuator  (model  presented  in  the  Task  4  Section)  is  used  to  modulate  the  bottom  mirror  of 
the  CCR. 

Figure  Tl-3  depicts  the  experimental  setup  for  determining  the  functionality  of  the  CCR.  Results  show  that 
the  angle  of  reflected  light  deviates  from  the  angle  of  incidence  by  less  than  one  degree  with  zero  applied  voltage. 
The  CCR  thermal  actuator  is  modulated  with  a  3  V  peak  square  wave. 
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Figure  Tl-2:  Captured  video  images,  from  different  perspectives,  of  solder  self-assembled  CCRs.  (Top)  200  x  200 
|tm  mirror  CCR.  (Bottom  Left)  500  x  500  p.m  mirror  CCR.  (Bottom  Right)  500x500  p,m  reflector. 


Microscope 


Figure  T 1-3:  Experimental  setup  for  veriiying  operation  of  the  CCR. 
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Task  2:  Array  of  3-D  Switches 


This  section  summarizes  development  of  solder-assembled  3D  switches.  Figure  T2-1  (a)  is  a  SEM  of  two, 
normally  open  (N.O.),  and  one,  normally  closed  (N.C.),  electrostatically  actuated  switches.  Figure  T2-1  (b)  is  a  top 
view  of  the  same  switches  showing  the  appropriate  operation  connections. 


Figure  T2-1:  (a)  SEM  of  a  normally  closed  (N.C.)  and  two  normally  open  (N.O.)  elctrostatically  actuated  switches  and 
(b)  a  top  view  of  the  same  switches  showing  electrical  connections. 

The  switches  were  assembled  with  4  mil  diameter,  63Sn/37Pb,  solder  spheres.  Each  plate  position  was 
controlled  using  mechanical  stops.  The  Signal-Ground  resistance  was  measured  at  6  ^2.  The  snap-through  voltage 
for  the  N.O.  switch  was  144  V,  and  55  V  for  the  N.C.  switch.  Each  switch  is  comprised  of  three  solder  assembled 
plates:  a  grounded  contact  plate,  a  signal  contact  plate,  and  an  attraction  plate.  The  N.O.  switch  is  operated  by 
applying  a  high  voltage  to  the  attraction  plate,  thereby  attracting  the  grounded  contact  plate  and  pulling  the 
grounded  contact  plate  against  the  signal  contact  plate.  The  N.C  switch  is  operated  by  applying  a  high  voltage  to 
the  attraction  plate,  thereby  attracting  the  grounded  contact  plate  and  pulling  the  grounded  contact  plate  away  from 
the  signal  contact  plate.  The  grounded  and  signal  contact  plates  are  held  in  contact  at  room  temperature  by  utilizing 
gold  on  polysilicon  residual  stress  in  the  grounded  contact  plate.  A  close  view  of  the  N.C.  switch  is  shown  in  Figure 
T2-2. 

The  solder  assembled  switch  designs,  based  on  three  vertical  solder-assembled  plates,  have  two  drawbacks. 
The  first  drawback  is  attributed  to  permanent/semi -permanent  charge  buildup  in  the  plates.  The  charge  buildup  may 
be  eliminated  in  the  future  by  using  alternative  materials  for  the  plates  and/or  by  proper  drive  electronics.  Secondly, 
the  optimum  combination  of  attraction  area  and  contact  plate  flexibility  is  very  important.  For  example,  some  designs 
have  relatively  low  attraction  voltages  because  the  area  of  attraction  is  large,  however,  the  grounded  contact  plate 
can  undesirably  come  in  contact  with  the  attraction  plate.  Optimization  of  attraction  area  and  plate  flexibility  is  no 
trivial  task  -  a  concentrated  effort  should  be  made  in  this  switch  design  if  better  results  are  desired. 
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Figure  T2-2:  A  close  view  of  the  N.C.  switch. 


Task  3:  Array  of  3-0  Microrobot  Less 

This  section  summarizes  the  development  of  3-D  microrobots,  A  solder  deposition  technique  using  an 
evaporator  has  been  developed  for  solder  assembly  of  microrobot  legs.  In  this  technique,  a  2  pm  thick  pure  indium 
layer  is  deposited  and  patterned  using  a  thick  photoresist  lift-off  process.  The  indium  is  reflowed  in  a  nitrogen  gas 
and  formic  acid  vapor.  Figure  T3-1  is  an  example  of  a  solder  self-assembled  microrobot  leg  array.  In  this  example, 
indium  was  deposited  and  patterned  into  118x118  pm  squares. 

Figure  T3-2  shows  a  single  solder  self-assembled  two-degree-of-freedom  microrobot  leg  with  close  views  of 
the  left  and  right  indium  solder  joints  and  the  electro -mechanical  linkage  isolating  the  right  solder  joint  from  the  leg. 
The  legs  did  not  assemble  to  90°  as  planned  in  this  design.  The  specific  reasons  for  not  achieving  90°  include  the 
deposition  of  more  indium  than  planned  due  to  the  imperfected  deposition/lift-off  process,  mechanical  resistance 
from  the  flexible  locking  arms,  and  the  extensive  unintentional  wetting  of  the  indium  to  the  gold  wiring.  These  sources 
of  error  are  correctable  and  can  be  eliminated  in  proper  design  and  batch  processing.  The  equivalent  of  a  sphere  of 
indium  with  diameter  of  37  pm  was  deposited  on  each  set  of  25x25  pm  gold  pads  for  the  microrobot  leg  shown  in 
Figure  T3-2.  The  outline  of  the  location  of  the  gold  pads  is  also  shown  in  Figure  T3-2.  The  design  planned  on  the 
indium  wetting  only  within  the  area  of  the  gold  pads,  however,  the  indium  jumped  its  boundaries,  wetting  extensively 
along  the  gold  wiring  attached  to  the  pads.  Again,  this  error  is  correctable  by  providing  breaks  in  the  gold  layer. 


Figure  T3-1:  A  solder  self-assembled  microrobot  leg  array 
after  indium  reflow. 
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Figure  T3-2:  A  solder  self-assembled  two-degree-of-freedom  microrobot  leg  showing  close  views  of  the  left  and  right 
indium  solder  joints  and  the  electro -mechanical  linkage  isolating  the  right  solder  joint  from  the  leg. 

Figure  T3-3  shows  a  single  solder  self-assembled  one-degree-of-freedom  microrobot  leg.  The  equivalent  of 
a  sphere  of  indium  with  diameter  of  15  pm  was  deposited  on  each  set  of  8x10  pm  gold  pads.  This  design  also  did  not 
lift  to  90°  for  the  same  reasons  stated  above  for  the  two-degree-of-freedom  microrobot  leg.  Figure  T3-3  also  reveals 
the  results  of  reflowing  small  amounts  of  indium  on  gold,  note  the  crystallized  indium  coating  the  gold  pads  and 
wiring  in  the  close  view  of  the  figure. 


Figure  T3-3:  A  solder  self-assembled  one-degree-of-freedom  microrobot  leg  showing  a  close  view  of  the  center 

indium  solder  joints  and  outlines  of  gold  pads. 
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Task  4:  Desisn  for  MEMS  Self-Assembly 


This  section  summarizes  work  on  the  designs  for  MEMS  self-assembly.  We  have  developed  MEMS  self- 
assembly  schemes  based  on  thermal  actuation,  stressed  cantilever  actuation,  and  electrostatic  actuation. 

We  have  developed  a  method  to  model  thermal  actuators.  First,  for  the  given  input  current,  the  temperature 
distribution  throughout  the  actuator  can  be  calculated  by  the  heat  transfer  model.  The  temperature  distribution  is 
then  used  as  an  input  to  a  finite  element  model  that  calculates  the  thermal  expansion  and  hence  the  bending  of  the 
actuator. 

For  exampfe,  the  average  temperature  (AT)  difference  between  a  hot  arm  and  the  cold  arm  of  the  vertically 
deflecting  actuator  depicted  in  Figure  T4-1  is  calculated  by: 

Ar  =  (AL,„,-AZ,,„„)/(«x4„,), 

where 

^,o.  =  loiT-T,)dx, 

a  is  the  coefficient  of  thermal  expansion,  and  T  is  temperature.  Model  results  show  that  the  average  temperature 
difference  between  the  hot  and  cold  arm  is  158°C.  Figure  T4-2  is  a  graph  of  the  simulated  temperature  distribution 
throughout  the  actuator’s  cold  arm  and  a  hot  arm.  Note,  only  half  of  the  actuator  is  modeled  due  to  the  symmetry  of 
the  design.  The  data  in  Figure  T4-2  are  then  used  as  input  to  a  finite  element  model  that  can  simulate  the  deformation 
of  the  actuator  as  shown  in  Figure  T4-3.  Simple  mechanisms  aBsembled  using  thermal  actuators  have  been 
demonstrated  with  the  comer  cube  reflector  and  microrobot  legs. 


Figure  T4-1 :  CAD  of  the  vertically  deflecting  thermal  Figure  T4-2:  The  calculated  temperature  distribution 

actuator.  The  scale  is  in  micrometers  and  corresponds  to  througout  the  vertically  deflecting  thermal  actuator’s  hot 
the  scale  in  Figure  T4-2.  and  cold  arm  when  supplied  by  2.5  mA  at  2.5  V. 


Figure  T4-3:  Simulated  deformation  of  the  vertically  deflecting  thermal  actuator  with  an  average  hot/cold  arm 
temperature  difference  of  158  °C.  The  predicted  tip  deflection  is  4  pm. 
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We  have  developed  analytic  and  finite  element  models  to  predict  the  deflection  of  the  tip  of  a  stressed  cantilever, 
such  as  the  cantilever  depicted  in  Figure  T4-4.  A  stressed  cantilever  is  a  polysilicon  cantilever  coated  with  a  0.5  jim 
layer  of  gold.  After  release  of  the  cantilever,  the  residual  material  stress  and  thermal  mismatch  of  the  gold  layer  on 
the  polysilicon  cause  the  cantilever  to  curl.  Figure  T4-5  shows  two  different  designs  that  use  stressed  cantilevers  to 
assemble  microrobot  legs. 


Figure  T4-4:  SEM  and  interferometric  microscope  measurement  of  a  stressed  cantilever. 


Figure  T4-5:  Microrobot  leg  assembly  with  stressed  cantilevers. 

Designs  using  stressed  cantilevers  for  self-assembly  of  microrobot  legs  have  been  successfully 
demonstrated  (Figure  T4-5).  Designs  using  stressed  cantilevers  for  assembly  of  a  CCR  have  also  been  implemented. 

Actuation  in  the  form  of  electrostatic  scratch  drives  (SDAs)  can  also  be  used  to  assemble  MEMS.  Figure 
T4-6  shows  one  such  design  that  uses  SDAs  to  pull  up  a  mirror.  Other  designs  include  a  mechanical  lock  to  keep  the 
MEMS  at  desired  position  after  being  pulled  into  place  by  the  SDAs. 
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Figure  T4-6:  SDA  assembly  of  a  mirror.  Left:  video  stills  of  a  single  mirror  being  assembled  using  an  array  of  scratch 
drive  actuators.  Top  right:  diagram  of  SDA  driven  assembly  design.  Bottom  right:  SEM  of  basic  SDA  assembly  for  a 

mirror. 


Task  5;  Design  for  Solder  Self-Assembly 

This  section  summarizes  the  progress  of  the  designs  for  solder  self-assembly.  A  model,  based  on  the 
surface  energy  minimization  of  liquids  {Surface  Evolver,  Brakke  K.A.,  Surface  Evolver  Manual,  Version  1.99,  June  1, 
1995,  University  of  Minnesota  Geometry  Center,  Minneapolis,  MN55455;,  has  been  completed  that  produces  a  very 
accurate  depiction  of  the  solder  behavior  [4-6].  Solder  type,  structure  size,  geometry,  and  material  properties  can  be 
input  and  the  model  will  calculate  the  angle  of  rotation  of  the  MEMS  structure.  Figure  T5-1  shows  a  completed 
solder  model  for  a  single  hinged  plate  rotated  to  45”  off  the  substrate. 


Figure  T5-1:  Solder  model  for  a  solder-assembled 
hinged  plate  rotated  to  45”, 

An  example  of  assembled  plate  is  shown  in  Figure  T5-2.  Figure  T5-3  illustrates  experimental  verification  of 
the  model.  The  error  bars  represent  the  sensitivity  of  the  model  to  the  variation  in  solder  volume.  The  (63%  Sn,  37% 


Figure  T5-2:  Example  of  a  solder-assembled  hinged 
plate. 
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Pb)  solder  balls  used  are  commercially  manufactured  from  0.004"  to  0.016",  +/-  0.001"  in  diameter.  This  is  a 
considerable  amount  of  variation  in  volume,  so  its  effect  is  important.  The  experiments  showed  that  the  model 
predictions  were  accurate  to  within  +/-  2  degrees. 

Figure  T5-2  is  an  example  of  a  structure  in  which  assembly  is  purely  controlled  by  the  solder  surface  energy 
minimization.  The  angle  of  rotation  can  be  controlled  fairly  precisely,  as  shown  in  Figure  T5-3;  however,  there  are 
undoubtedly  applications  that  would  require  more  precise  control.  Several  locking  mechanisms  have  been  developed 
and  demonstrated  in  which  the  angle  of  rotation  is  a  function  of  design  geometry.  Thus,  the  effect  of  solder  volume 
variation  is  eliminated,  and  the  accuracy  of  assembly  is  increased.  Two  locking  mechanisms  are  shown  in  Figure  T54 
and  Figure  T5-5. 

The  locking  mechanism  shown  in  Figure  T5-4  is  called  a  “slat  side-lock.”  The  central  plate  is  locked  in  place 
by  two  plates  that  trap  it  from  each  side.  Tabs  on  the  side  of  the  center  plate  fit  into  tapered  slots  in  the  trap  plates  for 
fine  alignment.  This  type  of  locking  structure  involves  a  sequential  set  of  plate  rotations  to  work  correctly.  During  the 
heating  process,  smaller  volumes  of  solder  will  become  molten  slightly  quicker  than  larger  volumes,  and  thus  can  be 
used  for  sequential  assembly.  Also,  different  types  of  solder  with  different  melting  temperatures  could  achieve  the 
same  effect. 

Another  way  of  mechanically  controlling  the  angle  of  plate  rotation  is  through  the  use  of  a  “kickstand  lock,” 
as  shown  in  Figure  T5-5.  As  the  solder  rotates  the  structure  out  of  the  plane  of  the  substrate,  a  rigid  bar  attached 
with  a  spring  rotates  with  the  plate  until  it  slides  into  a  hole  or  grove  in  the  plate,  thus  stopping  the  plate  movement. 
If  designed  correctly,  the  bar  would  not  only  stop  the  plate  rotation,  but  also  lock  it  in  position. 


Figure  T5-3:  Comparison  between  predicted  and  measured  final  angles  for  solder-assembled  single -hinged  plates. 
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Figure  T5-4:  Side-locked  plate.  Figure  T5-5:  Kickstand  locked  plate. 


A  study  was  done  to  gauge  the  effects  of  the  locking  mechanism  on  the  accuracy  of  MEMS  assembly.  It 
was  found  that,  depending  on  the  design  of  the  locking  mechanism,  the  repeatability  of  the  angle  of  rotation  varied 
from  +/-  0.3°  for  one  structure,  to  +/-  1.5°  for  another.  The  best  locking  mechanism  was  found  to  be  a  kickstand 
locking  mechanism  (Figure  T5-5)  with  mechanically  rigid  locks,  with  an  average  rotation  of  90.7°  +/-  0.3°  for  ten 
assembled  samples. 

A  second  modeling  method  developed  predicts  the  equilibrium  angle  of  a  MEMS  plate  by  approximating  the 
shape  and  volume  of  the  solder  joint  with  simple  geometric  shapes  as  shown  in  Figure  T5-6.  This  allows  the 
modeling  of  the  changing  solder  shape  by  an  equation  describing  the  torque  on  the  plate  from  the  surface  tension 
and  internal  pressure  of  the  solder,  and  an  equation  describing  the  changing  shape  of  the  solder  joint  for  any  angle. 
This  model  is  much  faster  than  the  Surface  Evolver  model  but  requires  training  to  be  as  accurate  [8].  Figure  T5-7 
shows  a  comparison  of  the  two  models  and  a  volume -normalized  plot  of  equilibrium  angle  and  solder  volume 
predicted  by  this  model,  useful  for  square  plates  of  an  arbitrary  size. 


Figure  T5-7:  Comparison  of  the  two  models  and  a  volume  normalized  plot  of  equilibrium  angle  and  solder  volume 

useful  for  square  plates  of  an  arbitrary  size. 


One  of  the  unknowns  of  the  solder  self-assembly  mechanism  is  what  are  its  limitations  when  reducing  the 
size  of  the  MEMS  structure.  This  question  was  examined  and  three  major  areas  of  limitations  were  identified: 
Fabrication  Limitations,  Modeling  Limitations,  and  Solder  Deposition/Metallurgical  Limitations. 

Fabrication  limitations  play  a  monumental  role  in  the  viability  of  reducing  the  aze  of  the  solder  self- 
assembled  mechanism.  For  example,  the  MUMPs  technology  used  in  our  research  has  minimum  fabrication  geometry 
of  two  microns  (http://www.memsrus.com/cronos/svcsmumps.html).  One  aspect  in  which  this  is  a  problem  is  in  the 
substrate  hinge.  The  play  in  the  hinge  that  results  from  the  minimum  geometry  limitation  allows  the  pin  of  the  hinge 
to  move  up  to  two  microns.  The  effects  of  the  hinge  play  on  different  size  structures,  as  well  as  locked  and  unlocked 
structures,  were  assessed.  The  data  show  that  for  the  larger  solder-assembled  structures,  the  hinge  play  effect 
results  in  a  lower  angle  error  than  in  the  smaller  structures.  Moreover,  the  mechanically  locked  assemblies  yield 
consistently  lower  error  in  the  equilibrium  angle. 

During  solder  assembly,  the  position  of  the  hinge  pin  is  most  likely  governed  by  some  physical  phenomena. 
If  this  phenomenon  can  be  understood  and  predicted,  then  future  designs  can  be  created  with  this  in  mind  to  further 
increase  alignment  accuracy  of  assembly.  For  this  reason,  extensive  modeling  and  experiments  were  performed  to 
understand  the  hinge  position  phenomenon.  Two  possible  factors  were  examined: 

•  Liquid  surface  energy  controlled  positioning.  Is  the  hinge  position  a  function  of  the  potential  energy  of  the 
system? 

•  Plate  structural  effects.  Does  the  plate  deformation,  solder,  and  residual  stress  effect  the  hinge  position? 

To  evaluate  the  surface  energy  effect  on  hinge  position,  Surface  Evolver  was  used  to  find  the  minimum 
energy  of  a  rotated  plate  system  in  which  the  hinge  position  is  varied  according  to  the  two-micron  fabrication 
limitations  of  the  MUMPs  process.  Figure  T5-8  shows  an  example  of  such  analysis  for  a  plate  designed  to  rotate  90" 
using  an  8-mil  diameter  63%  Sn  /  37  %  Pb  solder  sphere.  In  this  case,  there  is  a  definite  minimum  energy  position  for 
the  hinge  pin  {toward  the  substrate  and  toward  the  solder).  For  all  systems  modeled,  similar  result  was  consistently 
found. 
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Figure  T5-8:  (Left)  Surface  Evolver  predictions  of  surface  energy  for  the  range  of  possible  hinge  positions.  (Right) 
Diagram  of  hinge  positions.  The  lowest  energy  position  is  at  the  lower  right  position  (0.5,  -1.375). 


In  addition  to  surface  energy  effects  on  hinge  position,  there  are  many  other  fabrication  and  process  issues 
that  will  effect  the  angle  control  precision  such  as  solder  thermal  shrinkage  and  residual  stresses.  Figure  T5-9  shows 
the  effects  of  such  factors.  Both  hinges  in  the  figure  are  on  identical  structures  except  hinge  one  is  located  directly 
next  to  the  solder  pad,  and  hinge  two  is  farther  out  on  the  plate  (200  pm).  The  figure  shows  that  the  hinge  pin  in  hinge 
one  is  in  the  Surface  Evolver  predicted  position  but  that  for  hinge  two,  the  pin  has  been  pushed  to  a  slanted  middle 
position.  This  is  not  a  random  occurrence,  finite  element  (FE)  modeling  has  shown  that  for  this  case  the  residual 
stresses  in  the  polysilicon  have  bent  the  plate  backwards. 

To  evaluate  the  structural  effects  on  hinge  position  and  angle  precision,  FE  modeling  of  a  solder  assembled 
structure  has  been  performed.  The  primary  effects  analyzed  are: 


•  Solder  solidification  and  shrinkage  during  the  cool  down  step  of  the  reflow  process. 

•  Residual  stresses  in  the  structural  members. 

•  Hinge  placement  along  plate  edge. 

•  Kickstand  lock  placement  and  design. 

•  Plate  contact  with  substrate. 

•  Induced  stresses  from  inter-metallic  formation. 


Figure  T5-10  shows  one  such  result  for  a  450  x  400  pm  plate  rotated  to  90°  using  a  kickstand  mechanical  lock. 
The  structure  is  symmetric  (solder  is  in  the  center),  so  only  half  of  the  device  was  modeled,  and  the  solder  shape  was 
approximated  with  a  cube  for  simplicity.  In  this  case,  the  model  results  (right)  show  that  the  hinge  pin  will  move 
toward  the  kickstand  lock  during  deformation.  If  substrate  contact  is  added,  or  the  kickstand  lock  is  removed  the 
plate  will  warp  in  an  entirely  different  manner.  Figures  T5-1 1  shows  the  variety  of  warpage  phenomenon  that  can 
occur  in  the  same  plate  under  different  conditions  (in  this  case  substrate  contact  and  residual  doping  stresses).  All 
models  are  under  loading  due  to  thermal  stresses  resulting  from  thermal  expansion  mismatches  between  the  solder, 
gold,  and  poly  silicon.  The  warpage  results  are  displayed  in  a  contour  plot  of  the  deformation  in  the  out  of  plane 
direction.  In  addition  to  thermal  stresses  the  models  account  for:  (top  left)  residual  doping  stresses  and  substrate 
contact,  (top  right)  residual  doping  stresses  but  no  substrate  contact,  (bottom  left)  substrate  contact  but  no  residual 
doping  stresses,  (bottom  right)  no  residual  doping  stresses  and  no  substrate  contact.  From  this  study  we  see  that  the 
position  of  the  hinge  pin  and  warpage  of  the  plate  are  very  dependent  on  design  and  process  features.  Therefore, 
proper  design  should  minimize  or  eliminate  many  sources  of  plate  angle  variation  and  greatly  increase  the  overall 
assembly  precision. 
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Figure  T5-9;  SEM  photos  of  sample  hinge  Figure  T5-10:  Abaqus  FEM  of  deformation  occurring  due  to  solder 

pin  positions  that  can  result  from  varying  solidification  and  shrinkage.  (Left)  Defined  structural  areas.  (Right) 

design  parameters.  Resulting  deformation  scaled  by  a  factor  of  1 0. 
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Figure  T5-1 1:  Finite  element  model  of  a  solder  assembled  MEMS  plate.  Top  Left:  residual  doping  stresses  and 
substrate  contact.  Top  Right:  residual  doping  stresses  but  no  substrate  contact.  Bottom  Left:  substrate  contact  but 
no  residual  doping  stresses.  Bottom  Right:  no  residual  doping  stresses  and  no  substrate  contact. 


We  have  developed  and  experimentally  identified  FE  models  that  predict  the  plate  behavior  for  the  different 
process  factors  stated  above.  The  next  step  is  to  optimize  the  various  design  configurations  such  as  solder  pad 
shape,  hinge  position,  and  plate  geometry  to  yield  the  optimum  assembled  angle  precision.  With  the  many  variables 
that  can  be  altered  the  minimum  deformation  solution  is  very  non-obvious  and  potentially  unsolvable  analytically.  To 
resolve  this,  the  finite  element  models  discussed  above  were  integrated  with  optimization  algorithms  to  find  optimum 
designs. 

The  finite  element  software  ABAQUS  was  used  to  model  the  structure  and  extract  relevant  data,  and  the 
optimization  algorithm  NLPQL  was  used  to  optimize  the  variables.  The  optimization  process  is  as  follows:  starting 
with  an  initial  guess  for  the  variables,  the  model  results  and  the  derivative  of  those  results  with  respect  to  each 
variable  is  evaluated.  These  guesses  are  used  to  create  an  ABAQUS  input  file,  which  is  then  processed,  and  the 
results  extracted.  Since  this  is  not  an  analytical  function  to  be  optimized,  the  derivatives  are  found  using  a  simple 
finite  difference  method.  These  values  are  then  used  by  the  NLPQL  program  to  generate  a  new  prediction.  This 
process  is  repeated  until  the  solution  converges  to  some  desired  tolerance.  Since  this  process  could  run  for  many 
iterations,  it  was  important  that  an  efficient  model  was  used.  For  this  reason  Kirchoff  composite  shell  elements  were 
used  over  standard  20  node  3D  solid  elements.  The  resulting  computation  time  per  model  evaluation  was  reduced  by 
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approximately  95%.  The  average  time  to  converge  to  an  optimum  solution,  when  run  on  a  It's  500  MHz  DEC  Alpha 
(DEC  personal  workstation  model  500au)  with  768  MB  of  RAM  running  DEC  Unix  V4.0  was  about  30  minutes  and 
involved  ten  step  iterations  and  40  model  evaluations.  For  the  sake  of  brevity,  the  optimization  presented  here  is  for  a 
simple  plate  structure  and  for  a  basic,  building  block  component,  solder  self-assembly  structure. 

The  simple  plate  was  studied  to  gain  fundamental  insight  into  the  problem  of  plate  deformation  minimization, 
but  also  because  of  its  application  toward  designing  flatter  MEMS  mirrors  for  optical  applications.  Figure  T5-12 
illustrates  the  basic  structure,  a  two  microns  thick  polysilicon  plate,  point  fixed  at  the  center,  with  a  half  micron  thick 
gold  pad.  The  parameter  to  be  varied,  for  pad  and  plate  fixed  areas,  was  the  shape  of  the  gold  pad.  The  results  to  be 
optimized  were  chosen  to  be:  the  RMS  deflection,  average  deflection,  and  maximum  deflection  across  the  surface  of 
the  plate.  These  values  were  chosen  because  combined  they  best  depict  the  deformation  in  the  plate.  For  example,  it 
is  possible  that  a  plate  deformed  in  a  symmetric  way  could  have  an  average  deflection  of  zero,  but  by  no  means  be 
considered  flat.  Similarly,  a  plate  with  a  sharp  peak  in  a  concentrated  area,  could  still  have  a  low  average  and  nns 
deflection  value.  By  optimizing  all  three  values,  the  overall  parameter  to  be  optimized  is  indeed  the  flatness  of  the 
plate.  A  sample  result  is  shown  in  Figure  T5-13.  For  600  micron  by  300  micron  by  two  micron  thick  plate,  the  resulting 
optimum  pad  shape  resembled  a  t5rpical  dog  bone  shape.  This  result  is  very  non-obvious.  It  is  very  unlikely  that  the 
same  solution  could  have  been  found  with  a  standard  analytical  thin  film  plate  bending  model.  This  strongly 
illustrates  the  usefulness  of  using  finite  element  modeling  combined  with  and  optimization  algorithm  to  design 
MEMS  structures. 


Top  view 


Figure  T5-12:  (left)  Diagram  of  basic  plate  structure,  (right)  SEM  photo  of  basic  plate  structure. 


Figure  T5-13:  (left)  Sample  optimized  pad  shape  for  simple  plate  structure.  (Right)  Simulated  resultant  deformation 
where  red  is  ~1.5  microns  deflection  and  dark  blue  is  zero. 
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The  basic  building  block  solder  self-assembly  structure  consists  of  a  single  solder  sphere  with  a  hinge  and 
mechanical  lock  on  either  side.  A  typical  plate  is  between  100-500  microns  on  either  side.  A  600  micron  wide  by  200 
micron  tall  structure  was  chosen  to  optimize  because  ample  experimental  samples  were  available  to  confirm  the 
predictions.  For  this  case,  the  parameters  to  be  varied  were:  the  contact  position  of  the  mechanical  lock  and  plate,  the 
width  and  height  of  the  solder  pad,  and  the  position  of  the  hinge.  The  only  constraint  was  that  the  solder  pad  should 
remain  large  enough  to  be  practical  for  solder  deposition  processing.  The  plate  was  again  modeled  using  composite 
shell  elements,  but  the  solder  was  simulated  with  standard  three-dimensional  solid  elements.  The  interaction  of  the 
kickstand  and  hinges  was  modeled  using  contact  surface  approximations  rather  than  include  the  actual  hinge  and 
kickstand  structure  into  the  model.  This  greatly  increased  the  computational  efficiency  of  the  model  without  effecting 
accuracy. 

As  with  the  simple  mirror  structure,  the  optimization  program  was  able  to  generate  a  prediction  that 
significantly  optimized  the  deformation  in  the  plate.  Again,  the  values  to  be  minimized  were  the  RMS,  average,  and 
maximum  deflections  of  the  plate.  Figure  T5-14  shows  three  sample  cases:  a)  a  prediction  in  which  there  is  no  lock  or 
hinge  contact,  b)  a  prediction  in  which  the  lock  contact  position  and  hinge  have  been  placed  with  negative  results, 
and  c)  the  algorithm  prediction  for  lock,  hinge  position,  and  pad  dimensions  that  will  result  in  minimum  deformation. 
Interestingly,  the  case  in  which  the  lock  and  hinge  were  placed  poorly  resulted  in  a  more  severe  deformation  than 
with  no  lock  at  all.  The  poor  lock  and  hinge  position  resulted  in  a  maximum  deflection  of  ^5.5  micron  and  a  RMS 
deflection  of  ~3.4  microns,  whereas,  the  prediction  with  no  lock  or  hinge  resulted  in  a  max  deflection  of  ~4  microns 
and  an  RMS  deflection  of  ^2.1  microns.  Finally,  the  optimized  structure  showed  a  significant  improvement  with  a 
maximum  deflection  of --0.9  microns  and  RMS  deflection  of -0.6  microns.  If  the  max  deviation  angle  (resulting  from 
deformation)  were  calculated  as  the  deformation  from  flat  as  related  to  its  position,  the  average  angle  deviation  for  the 
optimum  structure  would  be  approximately  0.03  degrees. 


Figure  T5-14:  Comparison  of  results:  (top)  prediction  when  no  lock  or  hinges  are  used,  (middle)  prediction 
when  lock  and  hinge  position  are  poorly  chosen,  (bottom)  Algorithm  prediction  for  optimum  pad  size  and 

lock/hinge  position. 
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The  reason  for  the  reduced  deformation  is  likely  due  to  the  lock  and  hinge  constraints  working  against  the 
deformation  due  to  thermal  mismatch  with  the  solder.  The  solder  shrinkage  tends  to  cup  the  plate  around  it  like  a 
shroud.  By  placing  the  hinge  and  lock  near  the  edge  of  the  plate,  they  restrict  the  plate  and  force  it  back  toward  the 
desired  position.  If  the  lock  and  hinge  are  placed  too  close  to  the  solder,  they  only  amplify  the  deformation.  If  they 
are  too  far  out,  the  plate  will  bend  significantly  between  them  and  the  solder  joint.  It  should  also  be  noted  that  every 
structure  will  have  its  own  specific  optimum  design  that  depends  on  its  particular  geometry  and  processing.  It  is  also 
possible  that  for  some  applications,  global  flatness  may  not  be  the  optimization  goal.  For  example,  a  solder  assembled 
mirror  may  need  only  be  flat  in  a  specific  reflective  region.  Also  it  is  possible  to  maximize  curvature  for  added 
stiffness  in  one  direction  or  for  some  optical  application. 

The  above  work  discusses  the  reduction  of  plate  deformation  through  design  optimization.  However,  one  of 
the  primary  points  of  this  research  is  angle  precision.  At  this  point  it  should  be  reiterated  that  plate  deformation  and 
angle  precision  are  closely  related.  To  summarize,  because  the  hinge  position  of  both  the  lock  and  the  plate  is  known 
and  because  the  fabricated  dimensions  of  the  structures  are  very  concise,  copies  of  the  same  structure  should  all 
assemble  to  the  same  angle,  with  very  little  deviation.  But,  this  is  not  the  case.  Because  the  plate  deforms  due  to 
processing  stresses,  the  plate  will  deviate  from  the  intended  design.  The  question  is,  with  plate  deformation  how 
precise  can  angle  precision  still  be?  Using  a  mechanical  lock  it  has  been  demonstrated  that  0.5  degrees  of  variation  in 
angle  precision  was  possible. 

The  question  is,  using  the  optimized  design,  what  is  the  achievable  angle  precision.  Although  an  actual 
optimized  structure  has  not  been  tested,  similar  devices  with  slightly  less  than  optimal  designs  have  been.  In  one 
experiment,  six  400-by-400-by-3.5-micron  structures  were  assembled  and  measured.  The  pads  for  each  plate  were  130- 
by-130  microns,  and  the  hinge  and  lock  were  located  at  the  far  edges  of  the  structure.  After  measurement,  the  average 
angle  of  assembly  was  89.78.  More  importantly,  the  maximum  deviating  sample  was  at  89.69,  which  is  a  mere  0.09- 
degree  angle  deviation. 

The  position  of  more  complex  structures  made  from  multiple  solder  joints  has  also  been  investigated. 
Additive  rotation  and  position  error  occur  when  multiple  solder  joints  are  combined.  The  individual  solder  joints  and 
connecting  structures  (components)  have  characteristics  such  as  structure  dimension,  solder  volume,  temperature, 
hinge  play,  and  multiple  layer  induced  deformation.  All  of  these  characteristics  are  variable  and  may  or  may  not  be 
controllable. 

The  major  question  is:  given  a  set  of  individual  component  tolerances,  what  is  the  resulting  assembly 
tolerance?  This  question  will  become  increasingly  important  as  more  complex,  multi  part,  solder-assembled  devices, 
are  created.  For  example,  consider  the  three-hinge  optical  fiber  gripper  structure  shown  in  Figure  T5-15  (right).  The 
hinge  play  and  solder  variation  at  each  joint  will  combine  to  create  an  overall  position  error  at  the  end. 

To  investigate  the  overall  contribution  of  individual  component  tolerances  on  the  final  assembled  structure, 
a  model  has  been  created  that  takes  into  account  all  of  the  possible  variations  that  can  occur  in  a  complex  n -solder 
joint  solder  assembled  structure.  The  model  utilizes  standard  kinematic  analysis  combined  with  deformation  and 
solder  assembly  angle  models  to  find  the  assembled  shape.  The  range  of  error  for  the  final  assembly  is  determined  by 
applying  tolerance  analysis  techniques  to  the  model  (assembly  function  in  tolerance  analysis  terms). 

The  core  of  tolerance  analysis  is  to  find  the  assembly  function  that  relates  all  of  the  component  dimensions 
to  the  position  of  interest  on  the  assembly.  Figure  T5-15  (left)  is  a  sample  graphical  output  of  the  assembly  function 
for  the  fiber  optic  cable  gripper.  So  far,  the  model  requires  that  the  angle  variation,  hinge  play  effects,  structure 
geometry  variation,  and  temperature  (for  effect  of  bi-layer  residual  stress  induced  deformations)  be  inputted.  The 
purpose  of  this  tolerance  analysis  model  is  two -fold: 

1.  To  predict  the  tolerance  of  the  final  position  of  complex  solder  assembled  structures  given  existing 
process  capabilities. 

2.  To  investigate  the  effect  of  component  variations  on  the  assembled  structure. 
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Figure  T5-I5:  (Left)  Sample  graphical  output  of  the  assembly  ftinction  for  a  three-plate  fiber  gripper,  (Right):  SEM  of 
an  optical  fiber  gripper  assembled  with  8  mil  diameter  63Sn/37Pb  solder  spheres.  The  fiber  is  220  pm  in  diameter. 


This  work  is  essential  for  product  control  and  determination  of  design  feasibility  in  the  fliture  commercialization  of 
solder  assembled  microstructures.  Unfortunately,  a  model  is  only  an  estimation  of  the  effects  of  various  parameters 
on  assembly  control.  To  truly  understand  the  effect  of  various  parameters,  a  set  of  samples  were  fabricated  that 
would  study  the  actual  effect  of  structure  geometry,  internal  stresses,  pad  deformations,  hinge  play,  and  joint 
diameter.  The  test  structure  is  shown  in  Figure  T5-16.  The  samples  consist  of  seven  two  plate  assemblies,  each 
designed  to  assemble  to  different  angle.  For  example  the  first  structure  was  designed  such  that  the  first  solder  joint 
would  assemble  to  90  degrees  and  the  second  joint  would  assemble  to  45  degrees  and  is  thus  denoted  as  90/45.  The 
other  designed  variations  were  90/90,  90/135, 135/135,  135/90, 135/45,  and  45/135.  The  #-#  samples  of  each  device 
were  assembled  and  their  tip  positions  measured.  This  data  is  shown  in  figure  T5-17  plotted  against  the  model 
predictions  for  each  device  type.  It  was  found  that  solder  volume  variation  was  the  dominant  source  of  variation, 
followed  by  the  shape  of  the  solder  pad.  The  plate  shape,  internal  stresses  and  deformation  contributed  very  little, 
only  one  twentieth  of  the  impact  of  the  solder  volume  variation.  Interesting  the  hinge  play  contributed  the  smallest 
amount,  about  0.2  percent  of  the  overall  impact.  This  is  interesting  because  the  significant  amount  of  possible  hinge 
play  could  result  in  significant  variation.  The  reason  it  does  not  is  because  the  solder  consistently  pulls  the  hinge 
into  the  same  position. 


Figure  T5-I6:  Diagram  of  basic  test  structures  and  SEM  of  example  test  structures. 
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X  (microns) 

Resuhx  for  a  “90’, '4^’"’  axscmbty. 


X  (microns) 

Resulk  for  a  “90’, W’  asserobly. 


X  (microns) 

Resuhs  for  a  “OO’/nS’^’  assembly. 


Figure  T5-17:  Plots  of  sample  endpoint  measurements  versus  the  desired  position  for  the  seven  two-plate  angle 
combinations.  Lower  right:  sample  photo  measurement  of  a  typical  135/135  device. 


Solder  deposition  and  metallurgical  considerations  raise  the  following  questions:  how  will  the  solder  behave 
at  a  very  small  scale,  and  will  it  be  acceptable  for  MEMS  self-assembly?  A  test  was  performed  to  gauge  the  minimum 
size  in  which  surface  tension  still  dictates  the  shape  of  molten  solder,  barring  metallurgical  interactions.  A  1000- 
angstrom  layer  of  pure  indium  was  deposited  onto  a  non-wetting  silicon  wafer.  The  indium  was  then  reflowed  in 
formic  acid  gas  atmosphere  to  remove  and  prevent  oxidation.  Figure  T5-18  shows  that  nanometer  scale  solder 
spheres  were  very  well  formed.  In  order  to  see  the  smallest  of  the  solder  spheres,  it  was  necessary  to  use  an  atomic 
force  microscope  (AFM).  Using  the  AFM,  we  were  able  to  identify  solder  spheres  as  small  as  5  nm  in  diameter.  This 
gives  an  indication  that,  barring  metallurgical  interactions,  surface  tension  is  a  reliable  way  to  assemble  MEMS,  and 
that  the  surface  energy  model  is  still  valid  at  nano-scale. 
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Figure  T5-18:  Solder  nano-spheres  formed  by  surface  tension.  Left:  SEM  of  a  field  of  indium  micro  and  nano  scale 
spheres.  Middle:  Atomic  force  microscope  plot  of  a  field  of  nano-spheres. 


Other  considerations  in  MEMS  assembly  size  reduction  are:  can  the  solder  be  manufactured  smaller,  can  it 
be  deposited,  and  how  well  can  the  solder  volume  be  controlled?  These  questions  can  be  resolved  by  implementing 
a  thin  film  solder  deposition  system  for  deposition  directly  onto  the  MEMS. 

A  more  difficult  issue  is  metallurgical  limitations.  For  example,  molten  lead/tin  solder  will  absorb  gold,  which 
changes  the  material  properties  of  the  solder,  as  shown  in  Figure  T5-19.  For  larger  assemblies,  the  amount  of  gold  on 
the  pads  is  small  enough  to  not  change  the  material  properties  of  the  solder  significantly.  For  smaller  solder  volumes 
this  is  not  the  case,  the  solder  no  longer  behaves  in  an  easily  predictable  way,  as  shown  with  the  indium  on  gold 
pads  in  Figure  T3-3  and  Figure  T5-20.  Due  to  its  desirable  melting  temperature  and  wetting  properties,  the 
predominant  material  deposited  for  solder  self-assembly  is  indium.  Unfortunately  the  indium-gold  inter-metallic 
formation  rate  is  very  fast.  At  room  temperature  the  growth  rate  is  in  the  range  of  microns  per  hour.  This  means  that 
extensive  growth  has  occurred  even  before  the  devices  are  removed  from  the  deposition  chamber.  For  the  indium- 
gold  system,  the  vastly  predominant  inter-metallic  formed  is  Aulnz.  Auln2  is  brittle,  has  a  high  melting  temperature 
(457  °C),  and  has  material  thermal  expansion  properties  vastly  different  than  those  of  indium  or  gold. 

Although  some  research  has  been  done  on  the  inter-metallics  commonly  encountered  in  this  project,  there  is 
very  limited  data  available  on  the  diffusion  rates  and  inter-metallic  growth  rates.  For  that  reason,  we  studied  the 
existing  theory  by  experimentally  analyzing  the  inter-metallic  growth  patterns  in  the  MEMS  solder  self-assembly 
systems.  One  such  experiment  is  shown  in  Figure  T5-21.  Gold  disks  ranging  in  diameter  from  2  to  50  |im  were  coated 
with  an  equal  amount  of  indium  and  reflowed.  The  arrangements  with  a  gold  weight  percentage  larger  than  ~30  % 
showed  no  spherical  reflow  formation  (marked  50  and  45  in  Figure  T5-20),  indicating  complete  formation  of  Auln2. 
The  arrangements  with  a  gold  percentage  less  that  ~10  %  (marked  20  and  15  in  Figure  T5-20)  showed  complete 
spherical  formation. 


Figure  T5-I9:  Consumption  of  10  |im  gold  lines  by  an  8  Figure  T5-20:  Inter-metallic  effects  on  a  solder  sphere 


mil  Sn/Pb  solder  joint. 


approximately  25  |im  in  diameter.  The  indium  solder  was 


deposited  on  0.5  jim  of  gold  fabricated  through  MUMPs. 


22 


Figure  T5-21:  Picture  of  reflowed  indium  on  different  sized  gold  disks. 

Another  major  intermetallic  system  studied  was  the  lead-tin-gold  reaction.  This  particular  reaction  is  very 
well  understood,  as  tin -lead  (Pb-Sn)  solders  are  among  the  most  common  solders  used  in  industry.  In  the  tin -lead- 
gold  reaction  gold  is  readily  adsorbed  and  de-wetting  from  the  pad  may  occur.  While  the  room  temperature  diffusion 
may  have  long  term  effects,  it  should  not  greatly  effect  self-assembly.  Conversely,  the  molten  rate  of  consumption  is 
crucial  to  good  self-assembly.  If  all  of  the  gold  is  consumed  before  final  assembly,  the  solder  will  de-wet  and  the  final 
solder  shape  will  no  longer  be  what  was  intended.  Although  the  reaction  has  been  studied  thoroughly  in  industry, 
there  are  no  concrete  numbers  that  can  be  applied  to  solder  self-assembled  structures.  While  the  studies  speak  in 
generalities,  data  specific  to  our  particular  system  was  needed  to  gauge  and  potentially  control  the  de-wetting  effect 
when  using  Pb-Sn  solders.  To  do  this  a  simple  experiment  was  performed.  Four  and  eight  mil  63%  SN  27%Pb  solder 
spheres  were  placed  on  self-assembly  structure  pads  identical  to  those  used  for  actual  devices.  The  various  samples 
were  then  heated  and  analyzed. 

Figure  T5-22  shows  the  results  of  this  study.  From  both  sets  of  experiments  it  was  concluded  that,  when  the 
solder  wet  completely,  the  solder  would  remain  wet  for  at  least  the  first  1 7-25  seconds.  After  that  the  comers  of  the 
square  pads  begin  to  de-wet  and  pull  away  from  the  substrate.  After  two  minutes,  the  solder  has  completely  de-wet 
from  the  pads.  A  few  more  points  should  be  noted  about  this  data.  This  temperature  cycle  and  solder  were  chosen  to 
reflect  what  is  commonly  used  for  solder  self-assembly.  Other  temperature  ranges  and  solders  will  yield  significantly 
different  lesults.  Also,  the  length  of  time  each  solder  took  to  initially  wet  also  varied  (within  five  seconds  of  each 
other).  The  appropriate  measure  of  the  best  reflow  time  should  be  from  time  of  total  wetting  to  the  beginning  of 
de wetting.  Because  of  this,  these  experiments  have  some  degree  of  error  in  the  measurements.  For  this  reason,  it  is 
concluded  that  a  safe  maximum  time  to  allow  the  solder  to  stay  molten  and  reflow  should  be  less  than  ten  seconds. 
Unfortunately,  staying  under  ten  seconds  of  molten  time  is  not  easy.  In  addition  to  the  factor  stated  above,  the  rate 
at  which  the  stage  heats  and  cools  will  directly  effect  the  time  in  which  the  solder  remains  liquidus.  Controlling  the 
temperature  that  precisely  is  difficult.  To  effectively  use  tin-lead  solders  for  MEMS  self-assembly,  an  ultra-fast 


Figure  T5-22:  Left:  Total  dewetting  of  pad  (1  minute).  Right:  Total  wetting  of  pad  (10  seconds). 
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Intermetallic  reaction  is  an  extremely  complex  and  varied  problem.  Often  complicating  the  situation  is  the 
chemical  reaction  that  may  additionally  occur  due  to  the  necessary  use  of  flux  to  remove  metal  oxide.  An  example  is  a 
black  reactant  that  can  occur  in  the  area  surrounding  and  indium  solder  being  reflowed  in  formic  acid  vapor  flux. 
Under  the  right  conditions,  this  black  reactant  can  coat  adjacent  MEMS  structures.  This  reaction  is  a  particular  result 
of  the  combination  of  formic  acid  molecules  and  indium  oxide,  but  is  representative  of  a  very  common  set  of  problems 
that  arise  when  using  chemical  fluxes.  To  avoid  these  problems,  one  must  use  a  truly  fluxless  soldering  system.  We 
developed  a  prototype  of  such  a  system,  based  on  ultra -low  vacuum  removal  of  oxide.  The  system  is  based  on  the 
theory  that  if  the  partial  pressure  of  oxygen  is  low  enough,  the  Gibb’s  free  energy  of  the  oxidation  reaction  will 
change  such  that  the  oxygen  will  diffuse  off  of  the  metal  surface,  therefore  omitting  the  need  for  flux.  Figure  T5-23  is  a 
diagram  of  our  vacuum  soldering  system.  The  soldered  MEMS  are  placed  inside  a  temperature  controlled  vacuum 
chamber  that  is  itself  enclosed  in  a  chamber  filled  with  nitrogen  gas.  The  nitrogen  is  included  because,  at  the  very  low 
pressure  we  are  dealing  with,  it  is  near  impossible  to  prevent  some  leaking.  The  nitrogen  assures  that,  whatever 
leaking  does  occur,  it  will  be  nitrogen  that  leaks  in,  and  will  therefore  not  change  the  partial  pressure  of  the  oxygen. 
Successful  reflow  of  both  indium  and  tin  lead  solder  has  been  demonstrated  with  this  setup.  But  at  this  point  little  is 
known  about  the  varied  properties  of  the  solder  oxide  we  encounter.  Future  work  must  focus  on  quantifying  these 
properties  and  to  improve  the  system  for  MEMS  solder  self-assembly  in  a  vacuum. 


Figure  T5-23:  Diagram  and  photograph  of  fluxless  vacuum  soldering  system. 

Various  methods  currently  exist  in  industry  for  placing  solder  on  electronic  parts,  including  the  automated 
pick  and  placement  of  manufactured  solder  spheres,  molten  solder  droplet  printing  -  similar  to  inkjet  printing,  solder 
paste  printing,  evaporation  and  patterning  of  solder  metals,  and  electroplating  of  solder  metals.  In  our  research 
facilities  we  apply  solder  by  either  manual  pick  and  placement  of  individual  manufactured  solder  spheres  ranging  in 
size  from  2  to  16  mil  or  by  evaporation  and  photolithographic  patterning  of  pure  indium.  Solder  can  be  applied  and 
processed  on  the  MEMS  chips/wafers  before  or  after  oxide  release,  depending  on  which  is  more  convenient  for  the 
application.  Figure  T5-24  depicts  the  pick  and  placement  method  of  a  solder  sphere  and  an  example  of  patterned 
indium  before  reflow. 
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Figure  T5-24:  Depiction  of  the  pick  and  placement  method  of  a  solder  sphere  and  an  example  of  patterned  indium 

before  reflow. 

One  possible  way  to  eliminate  the  inter-metallic  formation  is  to  use  non-solder  materials  in  surface  tension 
based  self-assembly.  Figure  T5-25  shows  an  example  of  a  polysilicon  hinged  plate  that  has  been  assembled,  to  90° 
with  respect  to  the  substrate,  using  AZP4620  photoresist. 


Figure  T5-25:  Example  of  a  hinged  polysilicon  plate  assembled  to  90°  with  respect  to  the  substrate  using  AZP4620 

photoresist. 


AZP4620  is  a  positive  novolak  resin  based  photoresist.  The  resist  shown  in  Figure  T5-25  was  first  applied 
to  an  unreleased  MEMS  chip  and  then  patterned  and  developed  using  normal  processing  conditions,  where  soft 
baking  does  not  exceed  110  °C.  Subsequently,  the  chip  was  released  and  reheated  on  a  hotplate  in  air  up  to  160  °C 
whereupon  the  resist  shrinks  and  pulls  up  the  plate  over  a  period  of  approximately  5  minutes.  The  temperature  was 
reduced  once  the  free  plate  reached  an  angle  of  90°.  Any  positive  photoresist  should  exhibit  this  type  of  surface 
tension  and  shrinking  behavior.  For  AZP4620,  there  are  three  mechanisms  that  cause  the  shrinking,  and  these  are  as 
follows.  For  temperatures  above  120  °C,  the  remaining  PGMEA  solvent  evaporates  out  of  the  resist  and  the  DNQ 
begins  to  decompose,  causing  some  shrinkage.  Furthermore,  for  temperatures  above  125  °C,  the  three  diazo -groups, 
in  the  DNQ,  begin  crosslinking  with  the  novolak  resin,  resulting  in  the  tripling  of  the  molecular  weight,  causing 
shrinking.  Finally,  for  temperatures  above  150  °C,  the  novolak  resin  begins  crosslinking  with  itself,  also  contributing 
to  the  shrinkage.  Simultaneously,  surface  tension  forces  are  competing  with  the  curing  action  (shrinking)  causing  a 
relatively  slow  assembly. 

We  have  also  used  the  pure  shrinkage  of  200  pm  diameter,  uncrosslinked,  polystyrene  spheres  to  assemble 
plates.  When  heated  above  300  °C,  the  polystyrene  can  decompose  and  shrink  to  at  least  a  quarter  of  its  original 
size.  One  thing  to  keep  in  mind  is  that  the  polymers  that  we  have  experimented  with  are  not  appreciably  selective  as 
to  what  they  stick  to.  Our  observations  indicate  that  the  polymers  stick  equally  well  to  silicon,  silicon  dioxide,  silicon 
nitride,  and  gold.  Other  materials,  such  as  thiol-based  polymers,  which  attach  selectively  to  gold  should  be 
investigated  in  future  work. 
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Task  6:  Hybrid  Technology 


One  of  the  goals  of  this  research  was  to  combine  MEMS  actuation  modes,  with  solder  self-assembly  to 
investigate  solder  effects,  if  any.  We  have  combined  thermal  actuation  with  solder  assembly  in  the  solder  self- 
assembled  corner-cube  reflectors,  with  no  adverse  effects  noted. 

Another  hybrid  device  is  shown  in  Figure  T6-1.  A  released  surface  micromachined  silicon  chip,  containing 
an  array  of  six  resistive  point  heater  solder  self-assembly  devices  arranged  next  to  each  other,  was  packaged  and 
wire-bonded  to  a  14  pin  ceramic  dual  inline  package.  Volumes  of  63Sn/37Pb  solder,  each  equivalent  to  an  8mil 
diameter  solder  sphere,  were  deposited  across  the  gold  pads  of  the  hinged  plate  arrangement,  fabricated  over  the 
resistive  heater.  Power  (713  mW)  was  applied  to  the  resistive  heater  causing  the  solder  to  melt,  wet  the  gold  pads, 
and  assemble  the  surface  micromachined  plate. 


Figure  T6-1 ;  Solder  self-assembled  plate  heated  to  solder  melting  point  using  integrated  polysilicon  resistive  heater 

lines. 


We  have  also  combined  gold-on-polysilicon  bi-morph  cantilevers  with  solder  assembly  as  shown  in  Figure 
T6-2.  A  gold  covered  polysilicon  cantilever  is  used  to  connect  and  mechanically  isolate  the  solder  assembly  plates 
from  a  gold  mirror.  Because  of  the  nature  of  this  type  of  mechanical  isolation,  the  cantilever  was  covered  with  gold  in 
order  to  cause  the  cantilever  to  curl  upwards,  thereby,  pulling  the  two  connected  plates  up  a  little  to  initiate  assembly 
in  the  proper  direction  and  reducing  the  chance  of  mechanical  jamming  during  assembly. 

Finally,  the  solder  assembled  switch  is  an  example  of  electrostatic  actuation  combined  with  solder  assembly. 
We  have  also  demonstrated  an  electrostatic  rotary  scratch  drive  motor  with  solder  assembled  fins  attached  to  it 
(Figure  T7-9),  which  is  another  example  of  an  electrostatic/solder  assembly  hybrid. 
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Figure  T6-2;  Example  of  goId*on-poIysiIicon  bi-morph  cantilevers  combined  with  solder  assembly. 

Figure  T6-3  shows  a  hybrid  device  that  combines  micromirror  array  with  solder  assembly.  A  set  of  four 
electrostatic  micromirrors  is  lifted  parallel  to  the  substrate  to  achieve  greater  angle  of  rotation.  The  mirrors  rotate  +/* 
1 0  degrees  and  require  285  volts  to  operate.  The  entire  array  is  approximately  1  mm  across,  and  was  assembled  using 
8  mil  diameter  tin -lead  solder  spheres. 


Figure  T6-3:  Micromirror  array  combined  with  solder  assembly;  micromirrors  are  lifted  parallel  to  the  substrate  to 

achieve  greater  angle  of  rotation. 


Task  7:  Futuristic  MEMS 


This  section  summarizes  development  of  other  novel  MEMS  implemented  as  part  of  this  research.  Presented 
in  this  section  are  examples  of  multi-hinged  self-assembled  structures  that  have  been  demonstrated. 

Figure  T7-1  shows  a  polysilicon  box  enclosure.  The  box  is  525  pm  wide  by  525  |im  long  by  555  pm  high. 
Metalized  pads  are  125  by  125  pm.  8  mil  diameter,  60%In  /  40%Pb,  solder  spheres  were  used  to  assemble  the  box. 
Figure  T7-2  is  a  MEMS  garage  with  opening  door.  The  opening  door  is  attached,  at  its  base,  to  a  gear  that  is  driven 
by  a  thermal  actuator  array.  Figure  T7-3  is  a  14  hinged-plate  structure  that  was  assembled  with  three  solder  balls. 
Tethers  attached  to  the  tops  of  the  front  polysilicon  plates  are  used  to  lift  up  successive  rows  of  plates. 
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Figure  T7-1 :  SEM  of  a  polysilicon  Figure  T7-2:  A  MEMS  garage  with  Figure  T7-3  :  A  14  hinged-plate  structure, 

box.  opening  door. 

As  Figure  T7-3  shows,  the  concept  of  linking  MEMS  is  a  key  in  the  creation  of  more  complex  hinged 
structures.  Figures  T7-4  through  T7-6  demonstrate  this  concept  further.  Figure  T7-4  illustrates  two  types  of  mirrors 
assembled  using  solder  assembly  and  linkages.  Each  mirror  was  assembled  using  8  mil  diameter  63Sn/37Pb  solder 
sphere.  In  these  cases,  the  benefit  of  separating  the  mirror  from  the  solder  via  linkages  allows  for  movement  of  the 
mirror  section  via  external  actuation.  Figure  T7-5  is  a  photograph  of  five  plates  attached  end-to-end  and  assembled 
using  five  independent  8mil  diameter  solder  spheres  to  create  a  hook  or  arch.  The  device  was  designed  to 
demonstrate  the  complexity  that  is  achievable  with  linked  solder  self-assembly.  Each  plate  is  400  |xm  wide  and  450  pm 
long,  and  the  entire  structure  stands  approximately  750  pm  tall. 

Figure  T7-6  shows  a  design  of  a  linked  multi-plate  solder-assembly  fiber-optic  gripper/aligner.  The  device  is 
designed  such  that  when  the  solder  rotates  the  series  of  linked  plates,  the  motion  will  accurately  pull  and  lock  a 
single  optical  fiber  into  a  specified  position.  The  devices  in  Figure  T7-6  were  assembled  using  8  mil  diameter 
63Sn/37Pb  solder  spheres.  Figure  T7-7  is  a  design  of  a  solder  assembled  antenna.  The  device  is  100  pm  wide  by 
10,000  pm  high. 


Figure  T7-4:  Mirrors  assembled  using  solder  self-assembly  but  attached  via  linkages  to  allow  for  external  actuation 

methods. 
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Figure  T7-5:  Five-plate  hook  Figure  T7-6:  Linked-plate  solder-assembled  fiber  grippers  grasping  a  200  micron 
mechanism.  Each  plate  is  diameter  optical  fiber. 


400x450  |im 


Figure  T7-7:  Solder  assembled  antenna.  Device  stands  1  cm  tall  and  is  100  pm  wide  and  2jLim  thick.  It  was  assembled 
using  a  single  8-mil  diameter,  63%  Sn  /  37%  Pb,  solder  sphere. 

Through  the  combination  of  several  key  technologies,  we  have  demonstrated  a  solder  self-assembled, 
surface  micro -machined  axial  flow  fan  powered  by  a  scratch  drive  actuator  (SDA)  rotary  motor  (Figure  T7-8).  This 
micro  axial  flow  fan  has  the  potential  for  integration  into  chip  cooling  systems,  chemical  and  biological  gas  sensors, 
micro  combustion  mixing  chambers,  turbine  systems,  and  a  propulsion  system  for  a  micro -vehicle.  Not  only  is  the 
range  of  MEMS  applications  broadened  by  this  invention,  but  the  usefulness  of  some  existing  MEMS  components 
can  now  be  enhanced.  For  example,  rotor-stator  (wobble  or  toothed)  electrostatic  motors  can  be  outfitted  with  solder 
assembled  blades,  transforming  them  into  useful  fluid  propulsion  devices  instead  of  purely  laboratory 
demonstrations.  The  possibilities  continue  with  blades  mo  unted  on  comb  drives. 

The  axial  flow  fan  is  driven  using  a  rotary  SDA  array  made  from  nine  individual  SDA  with  a  continuous 
sliding  contact  around  the  perimeter.  Rotary  SDA  motors  have  been  previously  evaluated  in  the  literature  for  output 
torque  and  optimized  SDA  design,  yet  our  research  work  represents  the  first  demonstrated  application  of  such  a 
motor.  Each  of  the  eight  fan  blades  was  assembled  using  63Sn/37Pb  solder  with  volume  equivalent  to  a  4  mil  diameter 
sphere,  as  shown  in  Figure  T7-9. 

A  sample  size  of  27  blade  angles  was  measured  by  capturing  a  side  profile  of  each  blade  on  a  video  image.  The 
angles  were  calculated  from  the  pixel  distances  in  the  captured  video  images,  and  the  distribution  of  the  assembly 
angles  was  estimated  to  be  jiangie  “  107.24®  and  Oangie  =  7.95®.  The  SDA  design  and  drive  signal  have  been  optimized 
and  are  capable  of  driving  the  fan  at  speeds  over  180  RPM  in  an  ambient  air  environment.  The  frequency  of  the  drive 
signal  and  magnitude  of  the  driving  voltage  control  the  rotational  speed  of  the  fan,  as  well  as  the  direction  of 
operation.  Frequencies  from  1  to  12  kHz  were  used  to  drive  the  fan  with  a  50  V  zero-to-peak  sinusoidal  drive  signal 
offset  by  75  V  DC. 
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63Sn/37Pb  Solder  Ball 


SDA  Rotary  Motor 


Figure  T7-9:  The  top  shows  captured  video  images  of  an  unreleased  micro  axial  flow  fan  and  the  assembly  sequence, 


after  release,  while  solder  is  assembling  the  fan  blades  during  heating  on  a  hot  plate.  The  bottom  is  a  SEM  showing 


the  SDA  rotary  motor  and  solder  joints  with  an  exploded  view  of  the  fan  blade  hinges. 
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Task  8:  Quality  and  Reliability 


8.1  Introduction 

It  is  widely  recognized  that  warpage  occurs  in  multilayer  MEMS  structures  due  to  thermal  expansion 
mismatch  between  the  materials.  Figure  T8-1-1  (a)  shows  a  gold/polysilicon  plate  that  bends  about  20  pm  from  the 
substrate.  In  Figure  T8-1-1  (b),  a  gross  cylindrical  curvature  of  a  capacitive  plate  that  has  been  flip-chip  transferred 
to  another  substrate  is  seen. 


(a)  (b) 

Figure  T8-1-1:  (a)  SEM  photo  of  gold/polysilicon  plate  (MUMPS  28)  that  bends  upward  from  the  substrate,  (b) 
Capacitive  plate  (gold/poly silicon)  on  a  transferred  MEMS  structure  is  substantially  warped  due  to  residual  stresses. 

Figure  T8-1-2  shows  a  measurement  of  the  deformation  of  a  plate  after  it  has  been  assembled  (lifted  to  90°  off 
the  substrate)  using  an  8  mil  solder  ball.  The  image  is  taken  from  the  backside  of  the  assembled  plate  and  so  the 
solder  ball  is  hidden  from  view.  The  deformation  due  to  the  solder  self-assembly  force  is  evident.  In  this  case,  the 
observed  peak-to-valley  displacements  across  the  plate  are  3.7  pm. 


Figure  T8‘l-2:  Deformation  of  a  500  pm  x  420  pm,  2  pm  thick  polysilicon  plate  with  a  0.5  pm  thick  125x  125  pm  gold 
pad.  The  plate  was  assembled  with  an  8  mil  solder  ball  located  underneath  the  plate. 

This  section  summarizes  our  work  on  quality  and  reliability  of  developed  MEMS.  There  are  two  principal 
deformation  modes  that  we  studied:  i)  defonnation  of  multi-layer  plate  structures  due  to  thermal  expansion  mismatch 
between  individual  layer  materials  (e.g.  polysilicon  and  gold)  after  release,  but  before  solder  self-assembly;  and  ii) 
deformation  of  such  plates  after  they  have  been  assembled  using  solder.  The  experimental  efforts  consisted  of  fiill- 
field  measurements  of  deformation  over  the  entire  surface  of  plates  both  before  and  after  assembly.  These 
measurements  were  made  with  an  interferometric  microscope.  The  modeling  efforts  consisted  of  the  development  of 
idealized  models  of  multilayer  plate  structures  before  assembly,  that  consider  the  temperature  dependence  of  material 
properties  (stiffness  and  thermal  expansion).  Both  elastic  and  elastic -plastic  constitutive  behavior  is  included. 
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8.2  Experiment  capabilities 

We  measure  full  field  deformed  shapes  using  a  scanning  white  light  interferometric  microscope,  which  is 
shown  in  Figure  T8-2-1  (a).  Figure  T8-2-1  (b)  shows  the  optical  configuration  of  the  interferometric  microscope. 


(a)  (b) 

Figure  T8-2-1:  (a)  Scanning  white  light  interferometric  microscope  that  is  used  to  measure  the  full  field  deformed 
shapes;  (b)  optical  drawing  of  the  interferometric  microscope. 

A  custom-built  thermal  chamber  was  designed  that  can  be  used  together  with  the  interferometric 
microscope,  which  allows  to  measure  the  deformation  at  elevated  temperatures.  Figure  T8-2-2  shows  the  setup  used 
to  measure  deformation  at  elevated  temperatures. 


Ultra-low  load  indentation  can  measure  the  mechanical  properties  of  materials  at  very  small  scales.  The 
instrumented  microindentor  shown  in  Figure  T8-2-3  employs  high-resolution  sensors  and  actuators  that  can 
continuously  control  and  monitor  the  loads  and  displacements  on  an  indenter  as  it  is  driven  into  and  withdrawn  from 
a  material.  Many  material  properties  can  be  determined  by  analyses  of  the  indentation  load-displacement  data. 


Figure  T8-2-3:  Instrumented  microindentor. 
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8.3  Linear  and  geometrically  nonlinear  behavior  of  gold/polysilicon  fully  covered  plate  microstructures  subjected 
to  thermal  loading 

An  inherent  characteristic  of  multilayer  material  structures  is  that  misfit  strains  between  the  layers  (for 
example,  due  to  intrinsic  processing  stresses  or  thermal  expansion  mismatch  between  the  materials  upon  a 
temperature  change)  lead  to  stresses  in  the  layers  and  deformation  of  the  structure.  Curvature  measurements  are 
routinely  used,  e.g.  in  the  microelectronics  industry,  to  determine  stress  states  in  metal  films  deposited  on  a 
substrate.  Their  attractiveness  is  based  on  the  fact  that  wafer  curvature  can  be  easily  and  accurately  measured,  and 
through  use  of  the  Stoney  (1909)  formula  (based  on  small-deformation,  linear  elastic  considerations),  the  measured 
curvature  can  be  directly  related  to  the  film  stress  (which  is  typically  biaxial  and  spatially  uniform)  without  knowledge 
of  the  source  of  the  stress  or  even  the  thermoelastic  properties  of  the  film.  Much  of  the  understanding  regarding 
deformation  of  multilayers  that  has  come  from  microelectronics  applications  is  applicable  to  many  MEMS 
applications,  but  s^nificant  differences  exist  and  must  be  well  understood  to  optimize  the  design  of  reliable  MEMS. 
These  include: 

•  In  MEMS  applications,  the  thicknesses  of  the  layers  are  small  (on  the  order  of  |im)  and  usually  comparable. 
This  leads  to  much  larger  deflections,  relative  to  the  thickness  of  structures,  than  are  observed  in 
microelectronics  applications.  For  example,  the  maximum  deflection  of  a  0.5  jim  gold  film  on  a  500  pm  thick 
100  mm  diameter  sJicon  substrate  subjected  to  a  100”C  temperature  change  is  about  two  percent  of  the 
thickness.  These  values  are  not  unrealistic  for  microelectronics  applications.  For  MEMS  applications,  a 
reasonable  example  is  a  0.5  pm  gold  film  on  a  1.5  pm  thick  400  pm  diameter  polysilicon  plate.  For  this  case, 
when  subjected  to  a  100°C  temperature  change,  the  maximum  deflection  is  about  six  times  the  thickness. 

•  This  can  make  it  necessary  to  include  geometric  nonlinearity  in  order  to  accurately  model  deformation. 
Furthermore,  the  geometric  nonlinearity  can  lead  to  bifurcations  in  the  deformation  behavior.  These  can  be 
detrimental  when  dimensional  stability  is  a  requirement,  or  can  be  beneficial  for  actuator  applications. 

•  Since  the  layers  are  of  comparable  thicknesses,  stresses  can  vary  appreciably  through  the  thickness  of  the 
layers;  the  average  stress  in  the  layer  may  not  be  suitable  to  characterize  film  stresses  as  it  is  in  many 
microelectronics  applications. 

•  Tradeoffs  between  stress  and  curvature  exist;  for  a  given  metal  film  thickness,  decreasing  the  polysilicon 
thickness  can  reduce  the  stress  in  the  metal,  but  at  the  expense  of  increasing  the  curvature,  as  shown  in 
Figure  T8-3-1.  The  significance  of  this  is  obvious  as  many  MEMS  applications  have  strict  deformation 
requirements,  perhaps  more  severe  than  stress  requirements. 

•  The  curvature,  and  thus  stresses,  may  vary  significantly  over  the  in -plane  dimensions  of  the  structure;  thus, 
the  average  curvature  may  be  insufficient  to  adequately  describe  the  deformation  state  of  the  structure. 

•  The  patterning  geometry  of  metal  films  on  polysilicon  (or  other  materials)  may  be  arbitrary,  as  opposed  to 
fully  blanketed  or  thin  line  patterns.  This  can  result  in  complex  spatially  nonuniform  deformation  states. 


Gold 


Polysilicon 


Figure  T8-3-i:  (a)  Curvature  study  for  fixed  thickness  of  gold,  (b)  Stress  across  the  thickness  for  fixed  thickness  of 

gold. 

In  this  research  we  studied,  via  measurements  and  analysis,  the  deformation  behavior  of  a  series  of  square 
and  circular  gold/polysilicon  plate  microstructures  fabricated  by  the  MUMPs  surface  micromachining  process  and 
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subjected  to  uniform  temperature  changes  which  generate  internal  stress  and  deformation  via  thermal  expansion 
mismatch  of  the  gold  and  polysilicon.  Linear  and  geometrically  nonlinear  deformations  were  observed,  as  well 
bifurcations  in  the  equilibrium  deformed  shapes. 

We  designed  a  series  of  square  and  circular  gold/polysilicon  plate  microstructures  and  fabricated  them 
using  MUMPs.  The  square  and  circular  samples  were  implemented  using  the  MUMPS  3 1  and  36  fabrication  runs, 
respectively.  In  the  series  of  microstructures,  the  polysilicon  layer  was  fully  covered  by  the  gold  layer,  except  for  a 
couple  of  microns  around  the  perimeter.  We  varied  the  square  side  length  L  and  circle  diameter  D  to  include  L,D  = 
150  jj.m,  200  |im,  250  pm,  and  300  pm,  keeping  the  thickness  of  the  gold  and  polysilicon  fixed  at  nominal  values  of  0.5 
pm  and  1.5  pm,  respectively,  as  produced  by  the  MUMPs  process.  The  idea  behind  the  design  of  the  microstructures 
was  to  yield  a  gold/polysilicon  bilayer  structure  that  rests  as  freely  as  possible.  To  this  end,  the  plates  were 
supported  on  the  substrate  by  a  16  pm  diameterpolysilicon  post.  SEMs  of  typical  plate  microstructures  are  shown  in 
Fig.  Tg-3-2. 


Figure  T8-3-2:  SEM  of  L,  D  -  300  pm  gold/polysilicon  plate  structures.  The  support  post  and  etch  holes  are  apparent 
and  the  x-y  coordinate  system  used  in  subsequent  measurements  and  analysis  is  identified. 

We  measured  the  full-field  deformation  of  the  plate  structures  as  a  function  of  temperature  change  using  an 
interferometric  microscope  and  a  custom-built  thermal  chamber  that  is  covered  by  a  quartz  window  to  allow  optical 
access.  The  resolution  of  the  out-ofplane  displacement  measurement,  w(x,y),  is  0.1  nm;  the  resolution  of  the 
temperature  chamber  is  1°C.  The  specific  test  protocol  consisted  of  the  following: 

•  Heat  the  sample  to  approximately  100°C  where  it  is  flat  and  thus  presumably  stress  free.  The  sample  is  held 
at  this  temperature  for  a  time  period  sufficient  to  ensure  it  has  reached  thermal  equilibrium  and  a  stable  flat 
shape  (~  3.5  minutes). 

•  The  sample  is  then  slowly  cooled  to  room  temperature.  The  temperature  is  held  constant  for  ~  3.5  minutes 
roughly  every  5°C  so  that  thermal  equilibrium  is  reached.  Full-field  out-of-plane  displacements  yv(x,y)  are 
then  measured  over  the  entire  surface  of  the  plate.  This  data  is  stored  for  subsequent  analysis. 

•  From  the  measured  w(x,y),  we  calculated  the  curvatures  K^{x,y)~  -  d'^w/dx^  and  Ky{x,y)  ~  -  d^wjdy^  . 

These  are  computed  by  fitting  w(x.y)  with  a  6^-order  polynomial  in  x  and  y  (which  was  sufficient  to 
accurately  describe  the  displacement  profile  in  all  cases)  and  then  differentiating  this  function. 

For  each  sample  at  least  three  replicas  were  tested  and  excellent  repeatability  resulted. 

While  an  analytical  model  can  be  developed  to  model  the  deformation  of  two-layer  plate  structures,  the 
simplest  approach  to  tackle  more  general  deformation  problems  is  to  use  the  finite  element  method  to  solve  the 
geometrically  nonlinear  equations  over  an  arbitrary  spatial  domain.  This  is  also  the  most  viable  approach  for 
complicated  geometries.  We  used  the  ABAQUS  finite  element  code  with  composite  shell  elements  to  approximate  the 
thin-plate  kinematics  of  the  Kirchoff  theory.  Both  polysilicon  and  gold  were  modeled  as  linear  elastic  with  isotropic 
material  properties.  Input  parameters  to  the  finite  element  calculations  were  Ep=163  GPa,  Vp=0.22,  E^=78  GPa,  Vg=0.42. 
The  thermal  expansion  coefficients  of  the  materials  were  assumed  to  vary  linearly  with  temperature  and  values  at 
100(23)®C  used  were  otp=3.1(2.6)xlO'Vc,  and  otg-14.6  (14.2)xlO'V'’C.  Here  and  throughout  the  subscripts  p  and  ^ 
denote  polysilicon  and  gold,  respectively.  Calculations  were  carried  out  to  model  the  effect  of  the  support  post  and  it 
was  found  to  be  insignificant  on  the  resulting  displacements,  curvatures,  and  stresses  except  in  a  region  very  near 
the  post.  If  one  carries  out  calculations  as  just  described,  the  linear  and  geometrically  nonlinear  response  of  the  plate 
can  be  computed,  but  the  bifurcations,  and  the  subsequent  post-buckling  behavior,  can  not  because  of  the  perfect 
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symmetry  present  in  the  geometry  and  material  behavior.  In  order  to  model  the  bifurcations,  one  must  seed  an 
imperfection  into  the  model.  We  did  this  in  two  ways,  with  the  hope  of  identifying  the  likely  actual  imperfections  that 
contribute  to  the  observed  behavior.  First,  we  introduced  a  slight  thermal  expansion  mismatch  anisotropy  into  the 
model  by  assuming  the  thermal  expansion  coefficient  of  the  gold  is  orthotropic.  Specifically,  we  defined  the  in -plane 
thermal  expansion  coefficients  to  be  oc^  =  ttg  and  Oy  =  ttg+^a,  where  6a  was  taken  to  be  0.01%  of  Og.  With  this  slight 
perturbation,  the  pre -bifurcation  response  was  indistinguishable  from  that  with  the  isotropic  thermal  expansion.  The 
second  approach  to  seed  the  imperfection  was  only  used  for  the  square  plates;  we  defined  one  side  of  the  plate  to  be 
slightly  larger  than  the  other,  i.e.,  the  plate  has  dimensions  Lx(L+AL).  Details  of  the  bifurcation,  e.g.,  its  sharpness, 
are  strongly  influenced  by  small  changes  in  the  imperfection.  We  discuss  this  later  in  light  of  both  our  measurements 
and  predictions.  Calculations  were  carried  out  for  the  loading  situation  of  an  applied  uniform  temperature  change 
consistent  with  that  experienced  in  the  measurements. 

Figures  T8-3-3  and  T8-3-4  show  contour  plots  of  the  measured  and  predicted  displacement  w(x,y)  at  room 
temperature  for  all  square  and  circular  plates.  Each  contour  band  diagram  represents  a  displacement  of  0.23  |im, 
0.35|i.m,  0.45|Lim,  and  0.6  p-m  for  L  =  150  |im,  200  pm,  250fim,  and  300  pm;  and  a  displacement  of  0. 1 1 ,  0.24,  0.37,  and  0.43 
pm  for  the  D  ~  150,  200,  250,  and  300  pm  plates,  respectively.  Due  to  the  thermal  expansion  mismatch  between  the 
polysilicon  and  gold,  the  L,  D  =  150  jLim  sample  deforms  in  a  spherically  symmetric  manner;  contours  of  constant 
transverse  displacement  w(x,y)  are  nearly  circles.  This  is  also  the  case  as  the  size  increases  to  L,  D  =  200  pm  and  L,  D 
=  250  pm,  although  the  displacements  increase  as  the  plate  size  increases.  At  L,  D  =  300  pm,  though,  the  transverse 
displacement  contours  are  more  dliptical  than  circular,  indicating  that  the  deformation  is  no  longer  spherically 
symmetric.  It  is  apparent  that  upon  cooling  the  different  sized  plates  the  same  amount,  there  is  a  size  dependence  on 
not  only  the  magnitude  but  also  mode  of  deformation.  As  the  in -plane  dimension  of  the  plate  increases  with  the 
thickness  held  constant,  the  deformation  mode  changes  from  one  of  spherical  symmetry  to  one  more  like  cylindrical 
symmetry.  Both  the  measured  and  predicted  w(x,y)  contours  show  this  same  behavior  and  the  agreement  between 
them  is  quite  good.  As  we  will  describe  later,  proper  normalization  of  the  curvature  and  temperature  change  yields 
the  functional  form  of  this  size  dependence  before  bifurcation. 


(a) 


(b) 


Figure  T8-3-3:  Transverse  displacement  contours  at  room  temperature  for  square  gold/polysilicon  plate  structures  (L 
=  150, 200, 250,  and  300  pm  from  left  to  right):  (a)  measured,  and  (b)  predicted.  10  contours  for  each  contour  plot,  each 
contour  band  diagram  represents  a  displacement  of  0.23  p,m,  0.35p.m,  0.45pm,  and  0.6  pm  for  L  =  150  pm,  200  pm, 

250pm,  and  300  pm  respectively. 


35 


(b) 

Figure  T8-3-4:  Transverse  displacement  contours  at  room  temperature  for  circular  gold/polysilicon  plate  structures  (D 
=  150,  200, 250,  and  300  pm  from  left  to  right):  (a)  measured,  and  (b)  predicted.  10  contours  for  each  contour  plot,  each 
contour  band  diagram  represents  a  displacement  of  0. 1 1  pm,  0.24  pm,  0.37  pm,  and  0.43  pm  for  D  =  150  pm,  200  pm, 

250pm,  and  300  pm,  respectively. 

In  Figure  T8-3-5  we  plot  the  average  curvature  in  the  x-  and  y-directions  as  a  function  of  the  temperature  change 
during  cooling.  In  these  plots,  the  average  curvatures  are  determined  from  the  measured  and  computed  w(x,y)  by 

averaging  =- ^^w(x,0)/^x^  and  Ky  =-'d^w{0,y)/dy^  along  the  paths  y=0  and  x=0,  respectively,  over  a 

length  of  150  pm  from  the  center  of  the  plate.  The  x-  and  y-  directions  are  taken  to  be  aligned  with  the  principal 
curvatures  after  bifurcation.  Before  bifurcation,  the  response  is  axially  symmetric  and  so  the  x-and  y-directions  are 
arbitrary  and  indistinguishable.  The  use  of  the  average  curvature  as  a  measure  of  the  plate  deformation  is  obviously 
quite  appropriate  if  the  curvature  is  (or  is  close  to)  spatially  uniform.  This  aspect  will  be  taken  up  to  some  degree 
later.  The  measurements  and  predictions  in  Figure  T8-3-5  show  that  in  general  three  regimes  of  deformation  exist:  i) 
linear  thermoelastic  response  independent  of  plate  size,  ii)  geometrically  nonlinear  thermoelastic  response  that 
depends  on  plate  size,  and  iii)  bifurcations  in  the  curvature -temperature  response  that  also  depend  on  plate  size.  By 
bifurcations  we  mean  that  as  the  temperature  changes,  the  plate  develops  a  nearly  axially  symmetric  curvature  with 
~  which  at  first  increases  linearly  with  the  temperature  change,  and  then  nonlinearly.  At  some  critical  temperature 
change,  the  deformation  mode  changes  from  axially-symmetric  to  a  mode  closer  to  cylindrical  bending  where  k;.\ 
as  the  temperature  change  increases  beyond  this  value  the  difference  between  and  Ky  increases.  In  the  first  two 
regimes  the  measured  deformations  are  axially  symmetric,  while  in  the  third  they  are  strongly  nonsymmetric. 

Figure  T8-3-5  shows  good  agreement  between  the  measurements  and  predictions  in  all  three  deformation 
regimes.  The  major  discrepancy  is  the  sharpness  of  the  bifurcation  for  the  L,  £)  =  300  pm  plate;  it  is  quite  sharp  in  the 
predictions  but  much  more  gradual  in  the  measurements.  To  understand  this  we  remind  that  the  source  of  the 
bifurcation  is  an  imperfection  of  some  sort  that  breaks  the  ideal  symmetry.  In  the  calculations  that  imperfection  was 
simulated  via  the  slight  anisotropy  of  the  thermal  expansion  coefficient.  In  Figure  T8-3-6  we  explore  the  effect  of  the 
imperfection  that  triggers  the  bifurcation  on  the  deformation  behavior.  In  order  to  simplify  the  interpretation  of  the 
results  and  identify  basic  phenomenological  features,  the  calculations  in  Figure  T8-3-6  were  carried  out  with  constant 
room-temperature  material  properties.  This  slightly  alters  the  exact  magnitudes  of  calculated  quantities,  but  does  not 
effect  the  conclusions  drawn.  In  Figure  T8-3-6  the  thermal  expansion  coefficient  of  the  gold  ot^  =  ttg  and  Oy  =  a,+5cXj,, 
we  alter  the  imperfection  by  increasing  5ag  while  keeping  Wg  fixed.  Figure  T8-3-6  shows  that  increasing  5cXg  from  0.01% 
of  ttg  to  1.0  %  of  ttg  resulted  in  a  much  more  gradual  bifurcation  that  initiated  about  10  ‘’C  higher.  The  details  of  the 
bifurcation  are  strongly  influenced  by  slight  perturbations  of  the  imperfection.  This  sensitivity  will  probably 
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preclude  a  sharp  bifurcation  in  practice.  Indeed,  the  §aused  here  to  simulate  the  imperfection  is  not  an  unreasonable 
measure  of  the  anisotropy  that  can  develop  as  a  result  of  processing  variations.  It  is,  though,  probably  not  easy  to 
measure  the  thermal  expansion  anisotropy  this  accurately  with  available  means.  Furthermore,  other  imperfections 
such  as  slight  variations  in  the  circular  geometry,  the  etch  holes,  and  anisotropy  and/or  heterogeneity  in  other 
material  properties  also  can  play  a  role  in  the  exact  bifurcation  details.  As  a  result,  the  prediction  of  the  bifurcation 
details  is  prohibitively  difficult  because  of  the  strong  sensitivity  to  the  imperfection. 


100  90  80  70  60  50  40  30  20 
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Figure  T8-3-5:  Average  measured  (top)  and  predicted  (bottom)  curvature  as  a  function  of  temperature  upon  cooling 
from  100°C  to  room  temperature.  The  curves  from  top  to  bottom  are  for  the  L,  D  =  150, 200,  250,  and  300  pm 

structures,  respectively. 


Figure  T8-3-6:  Average  predicted  curvature  as  a  function  of  temperature  upon  cooling  from  100°C  to  room 
temperature.  L  =  300  mm  square  plate:  Sttg  from  0.01%  of  ttg  to  1 .0  %  of  ag. 

We  introduce  the  nondimensional  curvature  K  and  mismatch  strain  5: 
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Kf-  3AcATi^t,M. 

^  _ _ h  S 

In  these  expressions,  £  =  R  and  ^  =  L/2  for  the  circular  and  square  plates,  respectively;  is  the  layer  thickness, 
Mj  =  Ej  /(I  —  V^  )  for  both  the  polysilicon  and  gold  (i  =  p,g),  and  AoC—CCp  —  OC^  .  The  data  of  Figure  T8-3-4, 

before  bifurcation,  are  plotted  in  terms  of  these  nondimcnsionalized  variables  in  Figure  T8-3-7.  For  plate  structures 
with  in-plane  dimensions,  D  or  L,  much  larger  than  the  thickness,  the  average  normalized  curvature  should  be 
independent  of  plate  size  and  shape.  In  fact,  in  Fig  T8-3-7  predictions  are  shown  for  both  the  circular  and  square 
plates  and  they  are  indistinguishable.  The  measurements,  though,  differ  slightly,  particularly  at  larger  values  of  S. 
Here  the  results  for  the  circular  structures  are  slightly  above  those  for  the  square  structures.  The  reason  for  this  is 
unclear,  but  most  likely  is  due  to  slight  differences  in  properties  between  the  two  fabrication  runs  to  produce  the  two 
sets  of  samples. 


Nondimensional  mismatch  strain,  S 


Figure  T8-3-7:  Nondimensional  average  curvature  as  a  function  of  temperature  change  for  the  linear  and  nonlinear 
regimes.  Measurements  (open  symbols  for  square  structures  and  filled  symbols  for  circular  structures)  and  the 
prediction  (solid  line  for  square  and  circular  structures). 

In  Figure  T8-3-8  and  Figure  T8-3-9  we  explore  the  connection  between  the  thermomechanical  loading  (the 
temperature  change),  the  geometry  (plate  size),  and  the  boundaries  between  the  three  deformation  regimes  for  the 
gold/polysilicon  plate  microstructures.  Figure  T8-3-8  shows  the  temperature  change  necessary  to  initiate  nonlinear 
effects  (the  transition  between  regions  /  and  I!)  as  a  function  of  polysilicon  thickness  when  the  gold  film  thickness  is 
kept  constant  at  0.5  |im.  Figure  T8-3-9  shows  similar  results  for  the  onset  of  bifurcation  (the  transition  between 
regions  II  and  III),  For  simplicity,  the  calculations  leading  to  these  results  were  carried  out  using  an  analytical  model. 
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Figure  T8-3-8:  Temperature  change  required  for  the  Figure  T8-3”9:  Critical  temperature  change  for  bifurcation 
initiation  of  nonlinear  geometry  effects  as  a  function  of  as  a  function  of  plate  size  and  the  thickness  of 

plate  size  and  the  thickness  of  polysilicon  for  the  polysilicon  for  the  square  and  circular  gold/polysilicon 

gold/polysilicon  microstructures.  microstructures. 

The  FE  calculations  allow  the  q^atial  variation  of  the  curvature.  Our  fiill-field  measurement  allows  us  to 
study  this  experimentally.  In  Fig.  T8-3-10  the  curvatures  i^(x,y)  and  K(x,y)  are  plotted  as  a  function  of  temperature 
change  for  four  points  along  the  x-axis  (y==0)  for  the  L  =  300  pm  plate:  2x/L  =  0.17, 0.33, 0.50,  and  0.67.  Predictions  for 
these  four  locations,  plus  2x/L  =  1  are  also  shown  in  Fig.  T8*3-ll.  In  the  linear  regime,  the  curves  are 
indistinguishable  implying  that  the  curvature  is  uniform  across  the  plate.  In  the  nonlinear  regime,  though,  the 
curvature  varies  appreciably  with  position,  increasing  by  about  a  factor  of  two  from  the  center  to  the  periphery  of  the 
plate.  Good  agreement  exists  between  data  and  predictions,  with  the  predictions  slightly  overestimating  the 
curvature  in  the  nonlinear  regime.  The  spatial  variation  of  the  curvature  casts  doubt  on  the  suitableness  of  an 
analysis  based  on  constant  curvature.  The  most  widely  used  analysis  of  this  type  is  Stoney's  (1909),  which  is 
routinely  used  to  connect  measured  (uniform)  curvature  to  film  stress  in  film/substrate  systems. 
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Figure  T8-3-10:  Point  wise  curvature  as  function  of  temperature  change  for  the  L  =  300  pm  structure  at  four  locations, 
2x/L.  Solid  lines  are  predictions  and  the  open  circles  are  measurements. 
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Figure  T8-3-1 1 :  Predicted  pointwise  curvature  as  function  of  temperature  change  for  the  L  =  300  pm  structure  at 

various  positions  2x1  L  along  the  x-axis  (y=0). 

In  Figure  T8-3-12  we  plot  the  critical  average  curvature  at  which  the  bifurcation  in  the  equilibrium  shape 
occurs  as  a  function  of  the  plate  size  L.  Three  results  are  shown:  an  analytical  prediction  from  the  constant  curvature 
calculation  of  Salamon  and  Masters  (1995),  results  from  the  finite  element  calculations,  and  the  measurement  forthe  L 
=  300  |im  sample.  Despite  the  fact  that  the  constant  curvature  approximation  becomes  questionable  for  larger  plate 
sizes,  the  model  of  Salamon  and  Masters  for  the  critical  curvature  is  a  good  approximation  as  seen  by  the  good 
agreement  with  the  finite  element  calculations,  at  least  for  the  elastic  mismatch  and  plate  sizes  considered  here. 
Unfortunately  we  only  show  a  single  measurement,  but  it  is  accurately  described  by  both  the  analytical  and  finite 
element  results. 


Plate  Size  L  (pm) 


Figure  T8-3-12:  Critical  average  curvature  for  bifurcation  as  a  function  of  plate  size  for  the  gold/polysilicon 
microstructures.  The  solid  line  is  the  analytical  prediction  of  Salamon  and  Masters  (1995)  based  on  the  constant 
curvature  assumption.  The  open  symbols  are  predictions  from  the  finite  element  calculations  and  the  filled  symbol  is 

the  measurement. 


8.4  Linear  and  geometrically  nonlinear  behavior  of  gold/polysilicon  partially  covered  plate  microstructures 
subjected  to  thermal  loading 

Figure  T84-1  shows  measured  deformed  shapes  for  the  polysilicon  plates  with  square  gold  pads.  Without  a 
gold  pad,  there  is  a  small,  roughly  spherical,  deformation  due  to  processing-induced  intrinsic  stresses.  For  the 
smallest  gold  pad,  the  deformed  shape  is  also  roughly  spherical,  corresponding  to  the  shape  that  is  predicted  by  a 
linear  analysis.  However,  as  the  size  of  the  gold  pad  is  further  increased,  the  deformed  shape  changes  dramatically. 
For  the  two  larger  gold  pads  it  is  nearly  cylindrical.  Thus,  the  deformed  shape,  as  a  function  of  gold  pad  size,  exhibits 
a  bifurcation  as  the  pad  size  increases.  In  this  case,  the  deformed  shapes  exhibit  a  nearly  cylindrical  curvature. 
Figure  T8-4-2  shows  similar  measurement  results  for  the  polysilicon  plates  with  circular  gold  pads.  They  too  show  a 
bifurcation  from  spherical  to  cylindrical  curvature  as  the  pad  size  increases.  Interestingly,  though,  we  observe  what 
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appear  to  be  two  stable  modes  for  the  specimens  with  larger  gold  pad:  a  cylindrical  mode  deformed  across  the  flats  of 
the  polysilicon  plate,  and  a  cylindrical  mode  deformed  across  the  corners  of  the  polysilicon  plate.  These  deformed 
shapes  cannot  be  predicted  from  linear  theory. 


Figure  T8-4-1:  Measured  room  temperature  deformed  shapes  of  poly  silicon  (450  pmx  450  pm)/gold  (L  x  L)  test 
structures  as  a  function  of  gold  pad  dimensions:  a)  L  =  0  pm,  b)  L  =  100  pm,  c)  L  =  150  pm,  d)  L  =  200  pm 


(^)  (b)  (c) 

Figure  T8-4-2:  Measured  room  temperature  deformed  shapes  of  polysilicon  (450  pmx  450  pm)/gold  (D  =  diameter) 
test  structures  as  a  function  of  gold  pad  dimensions:  a)  D  =  0  pm,  b)  D  =  100  pm,  c)  D  =  150  pm 


In  order  to  better  understand  the  deformation  of  partially  covered  structures,  we  designed  a  series  of 
gold/polysilicon  partially  covered  plate  iricrostructures  and  fabricated  them  using  MUMPs.  A  CAD  layout  of  the 
test  structures  is  shown  in  Figure  T8-4-3.  All  test  structures  were  supported  by  16  pm  diameter  polysilicon  post  in  the 
center  of  the  plate. 


40  0  fim 

Figure  T8-4-3:  CAD  layout  of  test  structures  to  understand  and  guide  modeling  of  MEMS  deformation.  The  gray 
areas  are  poly  silicon  (1 .5  pm  thick)  and  the  teal  areas  are  gold  pads  (0.5  pm  thick). 

We  measured  the  full-field  deformation  of  the  partially  covered  gold/polysilicon  plate  structures  as  a 
function  of  temperature  change  using  an  interferometric  microscope  and  a  custom-built  thermal  chamber  that  is 
covered  by  a  quartz  window  to  allow  optical  access.  From  the  measured  y^(x,y),  we  calculated  the  curvatures. 
Measured  room  temperature  deformed  shapes  of  polysilicon  (300  pmx  300pm)/gold  (LxL)  test  structures  are  shown 
in  Figure  T8-4-4.  The  L=140  pm  structure  gives  a  symmetric  deformation.  The  L=160  pm  structure  shows  a  little  non- 
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symmetric  deformation  but  no  significant  bifurcation.  Bifurcation  occurs  for  L=I  80  )im,  L^200,  L=250,  and  L=300  jim 
structures. 


(d)  (e)  (f) 


Figure  T84-4;  Measured  room  temperature  deformed  shapes  of  polysilicon  (300  p.mx  300  pm)/gold  (Lx  L)  test 
structures  as  a  function  of  gold  pad  dimensions,  (a)  L  =  140  pm,  (b)  L  =  160  pm,  (c)  L  =  180  pm,  (d)  L  =  200  pm,  (e)  L  = 

250  pm,  (f)  L  =  300  pm 

A  finite  element  model  was  built  which  assumes  the  structures  are  stress  free  at  100  °C.  Measured  and 
predicted  results  showed  good  agreement.  Bifurcation  occurs  from  structures  with  L=180  pm  gold  pad  and  larger. 
The  curvature  development  in  x  and  y  direction  from  100  X  to  room  temperature  (25  °C)  are  shown  in  Figure  T8-4-5 
for  L  =  180  pm  and  250  pm.  Bifurcation  temperature,  the  temperature  at  which  the  curvatures  in  x  and  y  direction 
show  significant  difference,  increases  from  L=180  pm  to  L=250  pm.  But,  does  the  structure  with  a  larger  pad  always 
bifurcate  earlier  than  the  structure  with  a  smaller  pad?  The  answer  is  no.  Both  e^eriment  and  model  show  that  the 
L=300  pm  structure,  which  is  a  fully  covered  structure,  bifurcates  later.  Figure  T8-4-6  (a)  shows  the  bifurcation 
temperaUires  are  roughly  35  °C  for  L=180  pm  and  45  °C  for  L=250  pm,  respectively.  But  in  Figure  T8-4-6  (b),  the 
bifiircation  temperature  for  L=300pm  drops  to  about  35  °C. 
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Figure  T8-4-5:  Measured  and  predicted  curvature  development  during  the  cooling  process  from  1 00  °C  to  room 
temperature  for  the  center  supported  bilayer  structures.  The  bi  layer  structures  are  composed  of  polysilicon  plates 
(WxW,  W=300  pm)  and  partially  covering  gold  pad  (LxL):  (a)  L=180  pm,  (b)  L=250  pm. 
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Figure  T8-4-6:  Experiment  result  shows:  (a)  the  L=250  jLtm  structure  bifurcates  earlier  than  the  L=180  pm  structure,  (b) 
the  L=250  iLUn  structure  bifurcates  earlier  than  the  L=300  jLtm  structure. 


To  further  understand  the  curvature  development  of  partially  covered  structures  during  the  cooling  process, 
we  parametrically  studied  the  curvature  development  during  the  cooling  process  of  partially  covered  gold/polysilicon 
plate.  We  fix  the  size  of  polysilicon  plate  (300  |imx300  jjm)  and  vary  the  size  of  gold  pad  (Lx  L)  for  L=  100,  150, 175, 
200,  225,  250,  275,  and  300  pm.  The  predicted  results  are  shown  in  Figure  T8-4-7.  All  the  averaged  curvatures  are 
calculated  over  100  pm.  Figure  T8-4-7  (a)  shows  that  as  L  increases  from  100  pm  to  225  pm,  bifurcation  temperature 
increases.  Also  the  curvature  before  the  bifurcation  point  decreases.  But  as  L  increases  further  from  225  pm  to  300 
pm,  as  shown  in  Figure  T8-4-7  (b),  the  reverse  trend  occurs,  i.e.  the  bifurcation  temperature  decreases  and  the 
curvature  before  the  bifurcation  point  increases. 


(a)  (b) 

Figure  T8-4-7:  Prediction  results  of  the  cooling  process  from  100  °C  to  room  temperature  for  square  polysilicon  plates 
(300  pm  X  300  pm)  with  square  gold  pads  (LxL),  (a)  L=1 00,  i  50, 1 75, 200,  and  225  pm;  (b)  L-225, 250, 275,  and  300  pm. 

In  order  to  study  the  effect  of  the  peripheral  polysilicon,  we  compared  the  deformation  of  a  fully  covered 
structure  and  a  partially  covered  structure,  in  which  the  gold  pads  have  the  same  size  of  200  pm  x  200  pm,  but  the  size 
of  polysilicon  is  200  pm  x  200  pm  in  the  fully  covered  structure  and  300  pmx  300  pm  in  the  partially  covered  one.  We 
measured  the  average  curvatures  over  the  gold  pad.  Both  measurement  and  prediction  showed  that  the  fully  covered 
one  doesn’t  bifurcate,  but  with  a  peripheral  polysilicon  plate,  the  partially  covered  structure  bifurcated  for  a  -75°C 
temperature  change.  This  suggests  that  the  uncovered  circumference  polysilicon  can  significantly  affect  the 
deformation  behavior  of  the  center  gold/polysilicon  plate.  We  also  find  that  before  the  bifurcation  point,  the  partially 
covered  structure  has  lower  curvature  than  fully  covered  one.  This  suggests  that  the  peripheral  plate  restrict  the 
deformation  of  the  center  part. 


8.5  Thermal  cycle  experiment 
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Figure  T8-5-1  shows  the  curvature  variation  during  the  thermal  cycle  experiment  between  190  °C  and  room 
temperature  (24  °C)  of  a  center  supported  bilayer  gold  (0.5  \m\  thick)  /  polysilicon  (3.5  pm  thick)  square  plate  (L  x  L,  L 
=  200  |im).  The  heating/cooling  rate  is  about  10°C/3.5  minutes.  During  the  first  heating  cycle  thin  film  undergoes 
recrystallization  and  grain  growth.  There  is  no  further  grain  growth  in  the  subsequent  thermal  cycles.  As  shown  in 
Figure  T8-5-1,  the  curvature  development  repeats  in  the  temperature  cycles  after  the  first  heating  cycle. 


Temperature  (°C) 

Figure  T8-5-1:  Curvature  variation  during  the  thermal  cycle  process  between  190  °C  and  room  temperature  for  the 
center  supported  gold/poly  l+poly2  200  pm  x  200  pm  square  plate. 


8.6  Study  on  stress  relaxation 

For  many  applications,  it  is  important  to  accurately  control  the  deformation  of  multilayer  material  structures 
in  order  to  meet  the  performance  requirements,  especially  for  structures  subjected  to  thermal  loading  and/or  operated 
at  elevated  temperatures.  However  stress  relaxation  and  creep  can  significantly  influence  deformation  and 
compromise  structural  integrity. 

In  this  research  we  studied  both  experimentally  and  theoretically  the  stress  relaxation  behavior  of  a  series  of 
gold/polysilicon  MEMS  plate  microstructures  fabricated  by  the  MUMPS  surface  micromachining  process  and 
subjected  to  thermal  loading.  We  measured,  using  an  interferometric  microscope  with  a  custom-built  temperature 
chamber,  full  field  deformed  shapes  of  square  (I>:L)  gold  (0.5  pm  thick)/polysilicon  (1.5  or  3.5  pm  thick)  plate 
structures  of  sizes  L  =  200pm  and  L  =  300  pm.  The  structures  were  thermal  cycled  between  room  temperature  and 
190°C  to  stabilize  the  microstructure.  After  the  initial  thermal  cycle,  the  structures  were  cooled  from  190  °C  to  120  °C 
and  held  at  120°C  for  about  7  days.  We  measured  full-field  displacements  (and  from  these  determined  the  average 
curvatures  in  x  and  y  directions)  as  a  function  of  temperature  during  the  initial  thermal  cycle  and  the  subsequent 
temperature  drop  190  X  to  120  X,  and  as  a  function  of  time  at  120  °C.  The  deformation  during  the  temperature  drop 
results  from  thermoelastic  phenomena,  and  includes  linear  and  geometrically-nonlinear  behavior.  In  fact,  for 
microstructures  within  certain  size  and  thickness  regimes,  the  equilibrium  configurations  exhibit  a  bifurcation  from 
symmetric  to  strongly  nonsymmetric  deformation,  i.e.,  buckling.  During  the  isothermal  hold,  stress  relaxation  was 
observed  in  all  of  the  microstructures.  Interestingly,  for  microstructures  that  have  buckled  during  the  thermoelastic 
deformations,  the  stress  relaxation  can  cause  them  to  unbuckle.  Typical  results  are  shown  in  the  Figure  T8-6-1 . 
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Figure  T8-6-1:  Measured  curvature  development  during  cooling  and  a  subsequent  stress  relaxation  process  at 
elevated  temperature  of  a  300  |im  x  300  [im  gold  (0.5  ^m  thick)  /  polysilicon  (1 .5  |im  thick)  plate  structure:  (a)  Curvature 
development  during  cooling  from  190  °C  to  120  °C;  (b)  Curvature  development  due  to  the  stress  relaxation  when  the 

structure  was  held  at  a  constant  temperature  of  120  °C. 

We  attempted  to  model  the  stress  relaxation  process  by  assuming  simple  power-law  creep  ^  in  the 

gold,  and  assuming  the  polysilicon  did  not  relax  at  the  modest  temperature  of  120  °C.  We  found  that  a  stress 
component  of  n  =  7,  provided  self-consistent  results  among  all  measurements,  but  such  agreement  required  different 
values  of  A  for  the  different  thickness  microstructures  (Figure  T8-6-2).  The  microstructural  reasons  for  this  are 
unclear,  but  these  results  resonate  with  similar  studies  of  thin  films  on  thick  substrates  in  microelectronics  contexts. 
Interestingly,  in  our  samples  there  are  substantial  stress  variations  through  the  thicknesses  of  the  layers  since  they 
are  of  comparable  thickness. 


Figure  T8-6-2:  Measured  and  predicted  curvature  development  during  isothermal  stress  relaxation  process  at  elevated 
temperature  of  120  °C.  (a)  200  jimx  200  |a.m  gold  /  polysilicon  plate  structure  with  different  thickness  of  polysilicon; 

(b)  300  |amx  300  jum  gold/polysilicon  plate  structure. 
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8.7  Solder  deposition  test 

Indium  solder  was  deposited  on  top  of  the  gold  pad  of  the  gold/polysilicon  bilayer  structures  using 
evaporator  (Figure  T8-7-1).  In  order  to  study  the  intermetallic  reaction  between  indium  and  gold,  we  cast  the 
structures  into  epoxy  and  polished  the  cross  section.  Energy  dispersive  X-ray  (EDS/EDX)  analysis,  which  detects 
elements  by  their  emitted  X-ray  energy  signature,  was  carried  out  to  test  the  elements  on  the  cross  section.  Figure 
T8-7-1  (b)  shows  the  position  of  dot  test  and  line  scan  test. 


Figure  T8-7-1:  (a)  Multilayer  structures  after  indium  deposition;  (b)  schematic  diagram  of  the  plate  and  substrate,  red 
dots  indicate  the  locations  of  the  dot  test,  red  lines  indicate  the  locations  of  the  line  scan  test. 

The  line  scan  test  result  is  shown  in  Figure  T8-7-2.  Gold  and  indium  distributed  almost  at  the  same  area, 
which  indicate  the  intermetallic  reaction  between  gold  and  indium.  For  the  scan  length  of  6.076  pm,  the  intermetallic 
layer  is  about  1.5  pm.  The  dot  test  result  is  shown  in  Figure  T8-7-3.  Gold  and  indium  are  found  at  both  the  very  top 
of  the  plate  and  in  the  middle  of  the  plate.  One  dot  from  the  very  top  of  the  plate  shows  relative  equal  percent  of  gold 
and  indium,  which  suggest  compound  Auln  is  formed.  The  other  dot  from  the  middle  of  the  plate  shows  two  times  of 
indium  as  gold,  which  suggest  the  compound  Auln  2.  Dot  mapping  test  result  is  shown  in  Figure  T8-7-4.  This  again 
shows  intermetallic  reaction  between  gold  and  indium. 


(a)  (b) 

Figure  T8-7-2;  (a)  SEM  picture  of  the  line  scan  area,  (b)  Line  scan  result  which  shows  Si,  Au,  and  In  element 

component  distribution  along  the  scan  line. 
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Figure  T8-7-4:  Dot  map  results  which  show  element  distribution  at  the  area. 


8.8  Solder  reflow  test 

So  far  we  discussed  our  study  on  the  deformation  behavior  of  multilayer  structures  before  solder  self- 
assembly.  How  would  these  structures  deform  after  they  are  assembled?  We  did  the  solder  reflow  test  on  some 
partially  covered  structures  composed  of  circular  polysilicon  plate  and  circular  gold  pad,  as  shown  in  figure  T8-8-1 
(a).  The  solder  balls  are  60In/40Pb  alloy  and  the  diameter  is  5  mil.  During  the  reflow  process,  the  solder  will  spread 
only  on  the  gold  pad  and  the  shape  of  the  solder  after  reflow  is  formed  according  to  the  minimum  surface  energy 
principle.  At  about  205  °C,  the  solder  was  fully  reflowed,  and  then  the  structures  were  cooled  down  to  room 
temperature  in  the  reflow  chamber.  Figure  T  8-8-1  (b)  is  a  SEM  picture  of  room  temperature  solder  ball  after  reflow;  the 
gold  pad  size  is  200  pm  in  diameter  and  the  polysilicon  plate  is  450  jim  in  diameter.  Since  the  thickness  of  the 
polysilicon  plate  is  only  1 .5  pm,  it  is  very  thin  compared  with  the  dimension  of  the  center  solder. 
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(a)  (b) 

Figure  T8-8-I :  (a)  5  mil  solder  balls  were  put  on  top  of  circular  structures  and  re  flowed  ;  (b)  SEM  picture  of  room 

temperature  solder  shape  after  reflow. 


Although  the  center  solder  part  can’t  be  measured  under  the  microscope,  the  deformation  of  the  uncovered 
part  of  the  polysilicon  plate  can  still  be  measured  after  solder  reflow  test.  Figure  T8-8-2  (a)  shows  the  measured 
deformation  of  the  uncovered  part  of  polysilicon  plate  (diameter  D  =  450  jim)  after  the  reflow  process.  The 
circumference  polysilicon  plate  shows  a  3-fold  bifurcation.  This  happens  because  the  solder  increases  the  bending 
stiffness  of  the  center  part  and  thus  makes  the  bending  of  the  center  part  difficult.  But  as  the  temperature  change 
increases,  the  structure  can  not  deform  symmetrically  due  to  the  geometric  nonlinearity.  Thus  the  circumference 
polysilicon  plate  has  to  fold  around  and  the  S-fold  bifurcation  occurs.  Figure  T8-8-2  (b)  shows  the  predicted 
deformation  from  Abaqus,  which  quantitatively  agrees  with  the  experiment  result. 


Figure  T8-8-2:  Measured  and  predicted  deformation  of  polysilicon  plate  with  solder  in  the  center,  (a)  Measured;  (b) 

Predicted. 
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